ANALYTICITY OF THE PLANAR LIMIT OF A MATRIX MODEL 



STAVROS GAROUFALIDIS AND lONEL POPESCU 

Abstract. Using Chcbyshcv polynomials combined with some mild combinatorics, wc pro- 
vide an alternative approach to the analytical and formal planar limits of a random matrix 
model with a one-cut potential V. For potentials V{x) = — J2n>i o-nX^'/n, as a power 
series in all a„, the formal Taylor expansion of the analytic planar limit is exactly the formal 
planar limit. In the case V is analytic in infinitely many variables {a„}„>i (on the appro- 
priate spaces), the planar limit is also an analytic function in infinitely many variables and 
we give quantitative versions of where this is defined. 

Particularly useful in enumcrative combinatorics are the gradings of V, Vt{x) — a;^/2 — 
J2n>i^nt"^^x" /n and Vt{x) — x^/2 — J2n>3'^nt"'^'^~^x"' /n. The associated planar limits 
F{t) as functions of t count planar diagram sorted by the number of edges respectively 
faces. We point out a method of computing the asymptotic of the coefficients of F{t) using 
the combination of the wzb method and the resolution of singularities. This is illustrated 
in several computations revolving around the important extreme potential Vt(x) = x'^/2 + 
log(l — ^/tx) and its variants. This particular example gives a quantitative and sharp answer 
to a conjecture of t'Hooft's which states that if the potential is analytic, the planar limit is 
also analytic. 
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1. Introduction 

1.1. Formal Matrix Models and Their Planar Limit. Matrix models are integrals of 
exponentiated potential functions over finite dimensional vector spaces (such as the vector 
space of Hermitian matrices of size A^) that were studied in the seventies as an approximation 
of Quantum Field Theory in a 0-dimensional space-time. Matrix models at fixed value 
of N and their behavior when is large is useful in a variety of problems that include 
enumerative problems of ribbon graphs, random two-dimensional gravity, triangulations of 
surfaces, random matrices, topological string theory, intersection theory on the moduli space 
of curves and perturbative gauge theory; see [tH82, BIZ80, DFGZJ95, DV02, DV, Mar08, 
Mar, Meh04]. 

Matrix models come in two flavors: formal and analytic. Formal matrix models (FMM in 
short) are easy to define, using formal Gaussian integration. The input of a formal matrix 
model is a formal potential V 

(1) V{x) = ^-J2^x-eA[[x]], 

^ ft 

n=l 

which lies in a formal power series ring ^[[x]], where A is the completed ring 



(2) 



A = Q[[ai,a2,a3, ...]]. 
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The partition function Z and the free energy T of the formal matrix model is given by the 
following formal integral and its logarithm, respectively 

L (iMexp(-iVTr(V(M))) 

where 

• "Hat is the vector space of Hermitian matrices of size A^, 

• Tr(M) denotes the trace of a matrix M, 

• The meaning of the formal integration is the following: expand e-^Tr(v(M)+Af /2) 
formal power series in >1[[A^, Tr(M), Tr(M2 ),...]] and integrate coefficient -wise. This 
operation produces a well-defined element of iV^^[[l/A^^]]. 

So, we can write 



oo 

(4) J- = ^iV2-2^j-„ j-,eA 

5=0 



J-g G ^ is called the genus g-limit of the formal matrix model. We can expand J-g in terms 
of monomials 

A 

where the sum is over the set of all partitions X = (l'^i2"'2 • • • )' '^a = 11 j '^j^ ^^'^ '^a 
are rational numbers. J-g enumerates connected ribbon graphs of arbitrary valency on a 
connected, oriented surface of genus g; see [BIZ80, BIPZ78, LZ04, Pol05, LZ04]. More 
precisely, it follows by Wick's theorem that c\ is the weighted sum of all connected ribbon 
graphs (weighted by the inverse of the order of the automorphism group) of genus g that 
have Hk vertices of valency k; see [Pol05, LZ04]. When g = 0, J^q is the planar limit of the 
formal matrix model. The planar limit depends on the formal potential V, and if we want 
to stress this dependence, we will use the notation J^o^y. As an example, when 

2 

the coefficients of the formal power series J^o.v^ G Q[[a4]] counts the weighted sum of 
connected planar 4-valent ribbon graphs. From the definition of -Fo,V4) one can compute 
several terms of the power series J^oy^, by hand or by machine. The pioneering work of 
[BIZ80, BIPZ78] gave an exact formula for the power series J^o,Vi using potential theory: 

1 - 36a4 + 162a2 + (1 - 3004) VI - 12^4 1, f 1 - ^1 - 12aA ^„ 
^''^^ = 432^^ + 2 [ 6^^ ) ^ 

The computation of [BIZ80, BIPZ78] was lacking rigor, and several years later their method 
was justified by using potential theory and the Riemann-Hilbert method; see [EM03, DKM98]. 
In the present paper, we give an independent proof, using mostly techniques from real anal- 
ysis and elementary potential theory. In addition, we describe explicitly the analyticity 
properties of J-q with sharp results, see Theorems 1.1 and 1.2 below. 
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As a notational convention, we will use caligraphic symbols V,TZ,S,J^q,... for formal 
matrix models and straight symbols as V, R, S,Fq . . . for the analytic matrix models. 

1.2. Analytic Matrix Models and Their Planar Limit. Let us now define the analytic 
matrix models (AMM in short). An admissible potential V{x) is a function V : M — > M 
which is lower-semicontinuous, and grows sufficiently at infinity, i.e., satisfies 

(5) lim 4^ > 1- 

|a;|^-oo 2 log |x| 

For an analytic matrix model with an admissible potential V define 

(6) Iv = - lim A log f exp{-NTr{V{M)))dM 



inf < / V(x)fi(dx) — / / log \x — y\^(dx)^(dy) > 

-^vm [J J J J 



where V{M.) is the set of all probability measures on M. The second equality in the above 
equation follows for example from [Dei99], [Joh98]. 

In the case V{x) = ^ — Sn>i admissible potential we then define the analytic 

planar limit as 

(7) Foy = l-Iv. 

We will call Foy and ly the the analytic planar limit. 

As we already mentioned, this formula allows one to reduce the problem of the planar 
limit to the investigation of what is known in the literature as the logarithmic potential with 
external fields. 

A 1-cut potential V is an admissible potential whose equilibrium measure has support in 
a single interval [b — 2c, 6 + 2c]. Following the notation of [BDFG02], we will use the change 
of variables 

(8) {b,c') = {S,R). 



b + 2c b 6 -2c 5-2J?2 5 S + 2R^ 

Figure 1. A graphical interpretation of b, c and R, S in terms of the endpoints of the support 
of the equilibrium measure. 

It turns out that in the case V is a. 1-cut potential (plus some nondegeneracy) , and Va is 
an analytic perturbation of V, then the endpoints and the planar limit depend analytically 
on a. 

An admissible potential V is even if it satisfies V{x) = V{—x) for all x G M. For even 
1-cut potentials, the equilibrium measure of V is centered at 6 = 0. 
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1.3. Analyticity of the planar limit. Analyticity of functions in infinitely many variables 
is well defined and understood on functions defined on spaces (see [Lem99] and [Rya87]). 
In our case we need to define a weighted version of £^ space. To this end, let r > be a 
positive number, and for a complex-valued sequence a = {a„}„>i C C^, consider its £l norm 



(9) M\r = Yl 



n=l 



Now, consider the following type spaces 

= {a= {aj„>i C : ||a||,, < oo} 

4(2N) = {a = {aj„>i G iHN) : a^n = 0,n> 1}. 

Let Br and B^^ denote the open balls of radius 1 in and £j(2N), respectively. 

Now consider & C to be the set of sequences a = {a.„}„>i G such that 

2 n 



n 

n>l 



is a 1-cut admissible potential which is analytic near 0. Using Equation (7) we can define a 
map Foy 

(10) Fo : 6 — vR. 

Likewise, we have a map : — v R. 

We use here the definition of [Lem99, Rya87] for an analytic function on which 
essentially means that the Taylor series in infinitely many variables converges. 

Theorem 1.1. The maps Fq, R and S (resp. F^^ , and 5"^^j uniquely extend to analytic 
functions on -Bi/^ (resp. B^^^^). 

Our next theorem identifies the planar limit of the formal and analytic matrix model. 
Since the map Fq from (10) is analytic at G -^rl^)! its Taylor series regarded as a formal 
power series is given by 

(11) Fo = ^cxaxeA 

X 

where the sum is over the set of partitions (including the empty one), ca G Q and A is given 
in (2). Consider the formal power series {TZ,S) G A'^ defined in Section 1.4 below. 

Theorem 1.2. We have 

(12) R = n, S = S, Fo = To. 

What this means is that, if the analytical procedures are taken formally, one recaptures 
the planar limit of the formal matrix models. 

Theorem 1.1 and Theorem 1.2 confirm a conjecture of 't Hooft for the planar limit of 
matrix models, 't Hooft' s conjecture is motivated by perturbative gauge theory ideas whose 
Feynman diagrams are ribbon graphs, and asserts that J-o(V(x)) should be an analytic 
function at x = when V(x) is analytic at zero; [tH82]. For a proof of t' Hooft's conjecture 
for the case of Chern-Simons gauge theory, see [GL08]. 
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A natural problem is to extend Theorem 1.1 to all genera g. 

Problem 1.1. Show that for all g > 0, (resp., ^g^) is the Taylor series of an analytic 
function on B^/^ (resp. 5^/^). 

This may be achieved using [ACKM93, Eyn04]. 

1.4. Two gradings for the planar limit. The formal planar limit J^o & A enumerates 
planar ribbon graphs of arbitrary valency, and it is closely related to two other formal power 
series {TZ,S) which are uniquely determined by the system of non-linear equations 



n = ni(n,s) 



j qn-2j 



7j cn-2j-l 



(13) 
where 

(14) m{n,s) = 

n>l j>l J \ 3 ' 

(15) = E''"E("27)(j)''''' 

Equation (13) always has a unique solution in {TZ,S) G A"^ that satisfies i? G 1 + A~^ and 
5* G A'^, where A^ are the formal power series in the variables with no constant term. 
Moreover, it is easy to see that this unique formal solution has integer coefficients. 

An enumerative interpretation of the coefficients of {TZ,S) is given in [BDFG02], which 
in particular implies that they are natural numbers. An analytic interpretation of {TZ, S) is 
that they determine the position of the interval of a 1-cut analytic matrix model; see Section 
6. 

Enumerative combinatorics dictates two gradings on the set of variables a„, the edge 
grading degg(a„) and the face grading degj(a„,) defined by 

(16) dege(a„) = -, deg^(a„) = g ^ ^• 

Given an element G A, let "He G -^ll-^^^^]] and 7if denote the result of substituting a„ by 
anf^f"^ and ant"'^'^~^ respectively. For example, for the formal potential V{x) from Equation 
(1) we have 

2 oo 2 ^ 

(17) Ve(x) = ^r/V G A[[t'/\x]], Vfix) = G A[[t'/\x]] 

^ rt ^ rt 

n=l n=3 

where in the latter we assume that ai = 02 = 0. Likewise, for TZe, TZe and J^o,e- 
Of course, when we set t = 1 to Tie or "Hj, we recover Ti. In particular. 



(18) J-0,e(l)= J-0,/(l)= J-qG^ 

The next theorem gives a simple formula for J-o,e in terms of TZe and Se- This appears in 
[BIOS] but for polynomial potentials V and the proof in there uses orthogonal polynomials. 
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Theorem 1.3. We have: 

(19) = - f 'A^^mA^imum^,,, 

t Jq 2s 

It follows that 

(20) (enj = ^^Msm + ^'M, 

where f indicates the derivative with respect to t. 

The next theorem gives a simple formula for J-qj in terms oiTZf alone. 
Theorem 1.4. We have: 



(21) 



1 /■* 

^oj{t) = ^ {t-s)\og'JZf{s)ds 



In particular 

(22) it'Toj)" = \ognfit). 

Remark 1.2. Given a potential V = ^„>i an^^"/''^- as a formal power series, from the 
potential theoretic approach one obtains that at the formal level, c and b satisfy the system 



(23) 



n>l j>0 

In the case V = x^/2 — J2n>i o.nX"'/n, and R = c^, S = b, one easily obtains the system (13). 
Remark 1.3. Some authors prefer to consider the following rescaling Ve of Ve 



2t ^ n 

n>l 

Since Ve(a;) = Veit^^'^x), it is easy to see that 

c{t) = Vic{t) 
b{t) = Vtb{t), 

where c{t) and b{t) are defined using the system (23). 
Remark 1.4. Likewise, for the following rescaling V/ of V/ 
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we have Vf{x) = Vf{t^^'^x)/t which imphes that (here c{t) and b{t) are constructed from 
V{x)/t) 

c{t) = Vtc{t) 
hit) = Vib{t). 

Remark 1.5. When V is even, then S = and TZ satisfies the imphcit Equation 

(24) H{n) = l 
where 

^ oo 

(25) U{X) =X--J2 ('2nX 

n=3 

This is indeed so because 

1.5. Algebricity, holonomicity and asymptotics of the planar limit. In this section 
we discuss the algebricity of the planar limit. Let us recall first some well-known properties 
of algebraic functions and the asymptotics of their Taylor coefficients. The reader may 
consult [vdPS03] and also [FS09, Chpt.VII] for further details. Computer implementations 
are available at [DvHOl, Kau09, Pot07]. 

An algebraic function y = y{x) is one that satisfies a polynomial equation P[y,x) = 
for some 2-variable polynomial with rational coefficients. Below, we will be interested in 
algebraic functions y{x) which are regular at x = 0, i.e., they have a Taylor series expansion 

oo 

(26) y(x) = 5];a,x" 

n=0 

The set of algebraic functions is a field, closed under differentiation with respect to x. Al- 
gebraic functions are always holonomic i.e., they satisfy (regular singular) linear differential 
equations with coefficients polynomials in x with rational coefficients. An algebraic func- 
tion y{x) gives rise to a ramified d-sheeted covering C — > C with semisimple local mon- 
odromy (with eigenvalues complex roots of unity) and global monodromy a finite subgroup 
of SL((i, C). In other words, an algebraic function y{x) regular at x = can be uniquely 
analytically continued as a multivalued analytic function on C \ A, where A is a finite set of 
algebraic numbers. In practice the analytic continuation can be obtained via Puiseux series, 
and all local expansions of y{x) around a singularity x G A are exactly computed by ?/(x); 
see for example [DvHOl, Pot07]. Since y{x) is holonomic, it follows that the sequence (a„) of 
its Taylor coefficients from (26) is holonomic i.e., it satisfies a linear difference equation with 
coefficients polynomials in n with rational coefficients; see [Zei90]. To discuss the asymp- 
totics of (a„) we need to recall what is a sequence of Nilsson type, discussed in detail in 
[Garll]. 
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Definition 1.6. We say that a sequence (a„) is of Nilsson type if it has an asymptotic 
expansion of the form: 

(27) an ~„^oo >^"'n° (log nf Sx,a,l3hx,a,(3 

where 

• the summation is over a finite set, 

• the growth rates A are algebraic numbers of equal modulus, 

• the exponents a are rational and the nilpotency exponents (3 are natural numbers, 

• the Stokes constants S\^a,i3 are complex numbers, 

• the formal power series hx,a,i3{x) G are Gevrey-1 (i.e., the coefficient of is 
bounded by C"n! for some C > 0), 

• is a number field generated by the coefficients of hx,a,i3{x) for all A, a, p. 

For a detailed discussion of the uniqueness, existence and computation of the asymptotic 
expansion of a sequence (a„) of Nilsson type, see [Garll]. The results of [Garll] and the 
above discussion implies the following. 

Proposition 1.7. (a) If y{x) is algebraic and regular at x = 0, then the sequence (a^) 
defined by (26) is of Nilsson type, where /3 = in (27). 

(b) Moreover, the asymptotic expansion (27) can be computed exactly and effectively. 
We will apply the above proposition to the planar limit. 

Proposition 1.8. (a) If 7le(t),Se(t) (resp. lZf{t)) are algebraic functions, then if' J^'q e)' if) 
(resp. {t^ J^Qj)'" {t)) is also an algebraic function. 

(b) If V is a polynomial, then TZ^, and TZf are algebraic functions. 

(c) Let 

n>l 

Under the assumptions (a) it follows that the sequence (/„) is holonomic. 

(d) In addition, the sequence (/„) is of Nilsson type. 

Several illustrations of the above proposition to extreme potentials are given in Sections 
8 and 9. 

1.6. The plan of the paper. In Section 2 we introduce and discuss the formal matrix mod- 
els with the two important gradings, the edge and the face gradings. Section 3 introduces 
the potential theoretic part of analytic matrix models and the preliminary results needed 
in Section 4 where the main analytic results are presented. We use here real analysis tools 
combined with Chebyshev polynomials and elementary combinatorics to deal with the min- 
imization problem (7), which is an alternative to the classical complex analysis techniques. 

Section 4 is the bulk of the analysis, the central pieces being Theorems 4.2 and 4.3. These 
are applied to some analytic examples in Section 5. 

Next, in Section 6 we prove the matching claimed in Theorem 1.2 and in Section 7 we give 
the proofs of the main results for the formal matrix models, namely Theorems 1.3 and 1.4. 




10 



STAVROS GAROUFALIDIS AND lONEL POPESCU 



The main calculations with the extreme potentials are in Sections 8 and 9, for the edge 
grading and respective the face grading. These main calculations are complemented with 
a small discussion in Section 10 about the calculations in the case of planar diagrams with 
vertices of valence 3 or 4. 

In Section 11 we give the formal proof of t'Hooft's conjecture, materialized first in the 
general form of Theorem 11.1 and then in Corollary 11.1, from which Theorem 1.1 follows. 

At last. Section 12 gives a perturbation result which is used in the proof of Theorem 1.2 
in Section 6, though the results in this section do not give sharp results about the radius of 
convergence for the planar limit as in Section 11. However this is a very useful analytic tool 
and we decided to include here. 

Finally, the appendix contains some Taylor series of TZ, S and J-'. Some of these terms are 
used in the proof of the main results. Theorems 1.3 and 1.4. 

1.7. Acknowledgment. The authors wish to thank the anonymous referee for valuable 
comments, corrections and references. 

2. Formal matrix models 

It follows from the definition of the formal matrix model that the planar limit J-q is the 
generating series of counting of planar graphs, weighted by the inverse of the size of their 
automorphism groups. This is discussed in detail in [BIPZ78, Eyn, Mar, Pol05]. In particular, 
J-Q^e counts planar graphs where every ra-valent edge contributes a term t"/^. Likewise, J-Q^e 
counts planar graphs where every n-valent face contributes a term t"/^"^. 

3. Analytic matrix models 

3.1. A summary of analytic matrix models. One of the main problems one faces with 
the minimization problem (7) is the support of the equilibrium measure. Without extra 
assumptions on the potential V, the support can be an arbitrary compact subset of the 
reals. However, most of the formal computations on the planar limit as a counting object 
are based on the formal manipulations as if the support was one interval. 

Naturally, what we want to do here in the first place, is a complete analytical character- 
ization of the one interval support for the equilibrium measure of (7). The way we do this 
here is based on an elementary approach to the logarithmic potential due to the following 
formula for x,y E [—2, 2]: 

oo 2 

log|x-y| = -$:-T„(|)T„(|) 

n=l 

where T„ are the the Chebyshev polynomials of first kind. Based on this formula we give 
a quick incursion into various formulae in logarithmic potential theory on [—2, 2], especially 
the formula from Theorem 4.1 and show that the general case of one-cut potentials can 
always be reduced by rescaling and translation in the x-variable to this case. The reason of 
doing this is to highlight a way of using manipulations of the Chebyshev polynomials in this 
framework. This seems to be an alternative (in the case of measures with support [—2,2]) 
to the powerful complex analysis methods discussed for example in [ST97]. 
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However, the more interesting fact is that we obtain the following explicit formula for Jy. 
If y is a potential whose equilibrium measure has support [—2c + 6, 2c + 6] , then 
(28) 



, 2 vicx + h)dx r 

Iv = -\ogc+ , - / s 



2 TT\^— X^ Jo 



xV'{sx + b)dx\ ( f V'{sx + b)dx 



-2 27rV4 — x'^ J \J~2 vrv^4 — x 



ds. 



This is attained via concrete exploitations of the Chebyshev polynomials, first for the case 
of the interval [—2, 2] and then simple rescaling. It is worth pointing out that ultimately, 
this identity reduces to checking a combinatorial identity for binomial coefficients. This we 
carry out using the implementation of the Zeilberger method. 

The previous formula, makes the dependence on the potential very transparent. Any 
questions on the analyticity of /y (or F^y) under perturbation follows from the analyticity 
of the endpoints of the support of the equilibrium measure. 

Finally, we show that under certain non- degeneracy conditions made explicit in Section 12, 
if Vt is an analytic perturbation of V depending on the parameter t, then the planar limit 
Ivt depends analytically on t on a domain of the parameter space. 

Here is an outline of what follows. In Section 3.2 we introduce the main objects, in 
Section 3.3 we discuss the formula that connects the logarithmic potentials and the Cheby- 
shev polynomials. Next, in Section 3.4, we describe the connection with Fourier analysis. 
Section 4 contains the main analytical results. 

3.2. Logarithmic Potentials with External Fields. As it was pointed out in the In- 
troduction, we are going to look at the problem of minimizing the logarithmic energy with 
external fields and then investigate the planar limit in this framework. 

Assume ^ : M — )■ M is an admissible potential. For a closed set C M, according to [ST97] 
for the general case or [Dei99] for the case S* = M, the following minimization problem has 
a unique solution (which turns out to be compactly supported) 

(29) Iv{S) = M{Iv{fi):fieV{S)} 
where 7^(5*) stands for the set of probability measure on S and 

(30) /y(/i) = J V{x)fi{dx) - J J log \x - y\fi{dx)ii{dy). 

The term — // log \x — y\fi{dx)n{dy) is called the logarithmic energy of the measure /i. For 
simplicity, we will denote /y = /y(M). Also for a given measure /i, we will denote supp/i, 
the support of the measure fi. The equilibrium measure of (29) on the set S (cf. [ST97, 
Thm.I.1.3]) is characterized by 



V{x) > 2 / log \ x — y\fi{dy) + C quasi-every where on S 

(31) 

V{x) = 2 / log \x — y\fi{dy) + C quasi-everywhere on supp/x. 



Here, a property P holds "quasi everywhere" on the set Q if we can find a set A such that 
fi{A) = for any measure /i of finite logarithmic energy and the property P holds on Q\A. 
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This means, that the equahty on supp/x is almost surely realized with respect to any measure 
of finite logarithmic energy. 

Notice here that if we change the variable of integration to x — )■ cx + 6 and y ^ cy + b, 
where c 7^ 0, then, with 

/ic,f, = ((■ -6)/c)#/i, 

(■/c standing for the multiplication by 1/c), and for a given function 0, the push forward 
is the measure defined by (f)^fi{A) = : G A}) for any Borel measurable A. 

Therefore we have 

(32) = j V{cx + b)nb,c{dx) - j j log \cx - cy\nh^c{dx)nb,c{dy) 

= Iv{-c+b)-\og{c)ifJ'b,c) = Iv{-c+b)ifJ'b,c) — logC 

which in turn results with 

Iv = Iv{-c+b)-log{c) = Iv(-c+b) - log(c). 

3.3. Logarithmic Potentials and Chebyshev Polynomials. Recall that the Chebyshev 
polynomials of the first kind Tn{x) are defined by 

(33) T„(cos6') = cos(n6') 

see for example, [01v97]. Alternatively, they are given by the recursion relation 

Tn+iix) = 2xTn{x) - T„_i(x), To{x) = 1, Ti(x) = x. 

Tn are the orthogonal polynomials for the arcsine law l[_i^i](a:)^jy=p. The following lemma 
is due to Haagerup [Haa] and we reproduce the proof here for completeness. 

Lemma 3.1 (Haagerup). (a) For any real x,y & [—2,2], x ^ y,we have 

log|x-,| = -f:^T„(|)T.(f^ 

n=l 

where the series here is convergent on x ^ y. 
(b) If X > 2 and y G [-2, 2], we have 



log \x — y\ = log 



X 



+ Vx^^ 



n 

n=l 



where the series is absolutely convergent. 

(c) The logarithmic potential of a measure on [—2, 2] is given by 

(34) J log \x - ylfxidx) = - ^ j J T„ (^| j fx{dy) 

where this series makes sense pointwise. 

(d) The logaritmic energy of the measure /i is given by 

(35) JJ log \x - y\fi{dx)fi{dy) = ~Yl^(yJ ) /^(^^)) • 

In particular JJ log \x — y\fi{dx)fi{dy) is finite if and only if J2'^=i n (/ (f ) l^{dx)Y is finite. 
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Proof. We first point out tliat for any complex number z ^ 1, with \z\ = 1, one lias tliat 



(36) 



log(l -z) = -Yl 



n=l 



z'- 
n 



wliere we take tlie brancli of log on C\(— oo,0]. Now, write x = 2cosm and y = 2cosf, and 
observe 

„, \ , . f u + v\ . f V — u 

X — y = 2(cos u — cos v) = 4 sm | — z — I sm 



and hence, 



log \x — y\ = log 



2 sin 



u + V 



+ log 



2 sin 



V — u 



logll-e^^^+^^l + logll-e^^^-")] 

Re (log(l - e*("+'')) + log(l - e^^^'""))) 



n=l 

OO 



— (cos(n(M + v)) + cos(n(t; — u))) 



n=l 

OO 

E- 



cosinu] cos m; 



n=l 

OO 



n 



n=l 



For the case x > 2 and \y\ < 2, then write x = 2 cosh u = + e where u = log ^ 
and y = 2cosw, thus 

log \x-y\= log (e"(l - e-"+™)(l - e"""™)) 

= u + log(l - e-"+*^) + log(l - e""-™) 



+Vx'^-4 



U 



E 

n=l 



-e ""cosfrawl 



n 



For the second part, for given — 1 < r < 1 we introduce the kernel 



L,(x, := - 2] — T„ (^- j T„ 



n>l 
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This can be computed for x = 2 cos u and y = 2 cos v for u,v E [0, vr) with u ^ v a.s 

cosinu) cos(nv) 

n 

n>l 



— — (cos(?T,(u + v)) + cos(r2(u — v))) 
n 

n=l 

- —Re (e*'^("+^') + e*'"("-^')) 



n=l 



= (36) 
1 

" 2 



log|l -e*("+")| +log|l -e^("-^)| 

(log(l + — 2r cos(m + w)) + log(l + — 2r cos(m — v))) 



Next, for any 6, 



which results with 



4 > 1 + r - 2r cos > (2-2 cos 



1 + ^ 



1 + r 

log 4 > Lr{2 cosM, 2 cost;) > 2 log — ^ h log 1 2 cosm — 2 cos w| 



or for x,y E [—2, 2], x y, 



1 + r 

log4 > Lr{x,y) > 21og^ h log \x - y\. 

This combined with Fatou's lemma yields that 

^liin j Lr{x, y)fi{dy) = j log \x - y\n{dy). 

The rest follows. □ 

The first consequence of the above proposition is the computation of the well-known arcsine 
law of an interval; [ST97]. 

Corollary 3.2. If u{dx) = l[-^2,2]{x)^^^=^ is the arcsine law of the interval [—2,2], then 

'O, |x|<2 



log \x\>2. 



(37) j log \x - y\u}{dy) 

If yU is a signed measure on [—2, 2] with finite total variation and finite logarithmic energy, 
then 

(38) j log \x — y\ix{dy) = c almost everywhere for all x G [—2, 2] 



if and only if fi(dx) = 1 [-2,2] (x) ^^l^^^'^lt'^ ■ Here, almost everywhere is understood with respect 
to the Lebesgue measure. Additionally, the constant c must be 0. 
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Proof. Because the density of u is even, it suffices to prove (37) for x > 2 oi x E [—2,2]. 
Equation (37) follows from the lemma and the fact that the series in (34) is convergent and 
is convergent also in L^(l[-2,2](a^) ^^l_^2 )- 

For the second part, integrating (38) with respect to the arcsine law and exchanging the 
integration one obtains that c = 0. Now, using equality (34) we obtain that J (|) ^J^{dy) = 

for all n > 1 and thus fi and 1[_2,2] (x) have the same moments and consequently 

must be equal. □ 

3.4. A Connection with Fourier Analysis. Take the map 

e : [0,7r] 9 ^ ^ 2cos^ G [-2,2]. 
For a given measure fi on [—2, 2] we define /i = /i o 0, the measure on [0, tt] such that 

ji{A) = /i(e(A)) 

for every measurable set A in [0,7r]. The map fi ^ ft from measures on [—2,2] into the set 
of measures on [0, tt] is a one to one and onto. The advantage of using this comes from 

(39) an := j T„ (^|) fi{dx) = j cos{n9)il{d9), n > 0, 

which tells us that the "moments" of fi with respect to Chebyshev's polynomials are seen as 
the Fourier coefficients of a measure on [0,7r]. 

The following result is standard and we state it without proof. 

Proposition 3.3. Given a sequence {a.„}„>o, with ao = 1, the following are equivalent: 

(1) There exists a measure on [—2, 2] such that an = f T (|) fi{dx); 

(2) {an}n>o is a bounded sequence and 



2 



2 VT 



defines a nonnegative distribution, i.e. for any smooth nonnegative function (p : 
[-2,2] 

(^,0) >0. 

In particular, if Yl^=i\'^n\ is convergent, then there is a measure fi on [—2,2] such that 
an = J T (l) fi{dx), if and only if 

oo 

(40) m(x) := l + 2^a„T„ (I) > for all xG[-2,2]. 



n=l 



In this case the 



measure fi is given by jj{dx) = -^^^ 



u{x)dx 
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4. The planar limit of analytic matrix models, The Main Results 
Given a continuous function / on [—2,2], we define 




n > 0. 



Notice that if / is bounded by a C^~^ function and piecewise C'^ for some k > 0, using 
the Fourier interpretation and repeated integrations by parts we learn that (3n{f) = o{n~'^). 
Next define the orthogonal polynomials 



(42) r„(x) = V2Tn{x/2) 

for n>l and Tq = Tq = 1. These provide a Hilbert basis of ^l[_2.2](x) ^/j-^ '^^^ 

Theorem 4.1. Assume that V is a C"^ and piecewise function on [—2,2] and A G M a 
constant. Then, there is a unique signed measure fi on [—2,2] of finite total variation which 
solves 

2 J log \x — y\fi{dx) = V{x) + C almost everywhere for x G [—2, 2], 
_/i([-2,2]) = A 

where almost everywhere is with respect to the Lebesgue measure on [—2,2]. The solution fj, 

T / I \ u(x)dx 

IS given by n[ax) = -^^^ — 



where 



m u{x^ = A-- r y^'^y^^y n^M^^i^ f v\x) - v\y) dy 

2y_2 7rv/4^ 2y_2 7rv/4^ J-2 a; - y ti^A^' 

In addition, the constant C must be given by C = — j'^^ 

Proof. In the first place, the uniqueness is clear because of Corollary 3.2. 
To prove the rest we first write the function V 

oo oo 

V{x) = Y,{f^,V)Ux) = (5o{V)+2j2(^n{V)Tr, (I) 

n=0 n=l 

where the inner product is taken in ^l[_2,2](a^) point out that the regularity 
of V implies that Pn{V) = 0(l/n^). Invoking representation (34), results with 

-2 E ^ ( / ^« (f ) /^(^^)) (f ) = ^ + ^o(^) + 2 E (^r.iV)T^ 

n>l ^"^ ^ n=l 

Thus, equating the coefficients, we must have now C = —f3Q{V) and 



J T„ (I) /i(rfx) = ~^n{V) 
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which, means that n{dx) = for 

oo 

u{x) = A-J2nPn{V)T^{£). 

n=l 

Here is the point where we need the assumption to make sure this series converges 
absolutely since in this case nf3n(y) = o(l/n^). 

To prove equality (43), our task now is to show that 

- Vn/3 (V)T (-] = -- C VYIM^-- C + C - V'{x) dy 



n=l 



Notice that both sides of this equation are linear functions of V and thus by a simple 
approximation argument it suffices to check it for the case of V{x) = (|) for some 
m > 1. After making the change of variables x = 2cos'U, y = 2cosf, this identity reduces 
to checking that 

1 /"'^ cosf sinnu cosu f"^ sinnv sinu s'minu) s'mv — sininv) s'mu 

— cos{nu) = / dv / dv -\ / dv. 

vr Jo sin v jo sin v ^ Jo (cos u — cos v) sin v 

Now, instead of checking this we look at the generating functions of the right and left hand 
sides. Specifically, using the fact that for — 1 < r < 1 and t G [0, vr], 

oo oo . 

, , , , « 1 /X COS t — r , r> 1 . / X sm r 

(44) > r" ^ cos(?it) = and > r" ^ sm(nt) = -, 

1 — 2r cos t + r-^ ^-^ 1 — 2r cos t + 

n=l n=l 

the rest follows from straightforward computations. □ 

Proposition 4.1. If fj, G P([— 2,2]) and V is a and piecewise function on [—2,2], 
then 

(45) Jv(/i) = (3o{V) + 2^^ bn{V)ar, + ^ j 

n=l ^ ' 

where 

(46) an = j Tn fi{dx). 
Furthermore, we also have that 

^ oo 

(47) Ivifi) > f3oiV) - -J2nf3n{V)' 

^ n=l 

with equality if and only if an = —/3n{V)/2 and 

oo 

(48) l-^n/3„(\/)T„(|) >0 for any xG[-2,2]. 

n=l 

In this case, 

dx 



(49) ^i(dx)= l-5^n/3„(l^)T„,(| 



n=l 



7r-\/4 — 
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Proof. One can write 

/oo « oo 

Vdfi = (3o{V) + 2j2Pn{V) / T„ (I) i2{dx) = (3o + 2 
n=l n=l 

To prove (47), one needs to complete the square in (45) to get that 



oo CO , 

iy[^) = /3o(y) - - ^ n^yf + E - ( 

n=l n=l ^ 



This imphes inequahty (47). The equahty is attained only for the case a„ = —Hih^^Yl which, 
cf. (40) determines a measure on [—2, 2] if and only if (48) is satisfied. The rest follows 
easily. □ 

We arrive at a necessary and sufficient condition for deciding that an equilibrium measure 
on [—2, 2] has full support. 

Corollary 4.2. Assume that \^ is a and piecewise function on [—2,2]. Then, the 
equilibrium measure on [—2, 2] has full support if and only if 

oo 

(50) 1 -^n/3n{V)Tn > for x on a dense subset of [-2, 2]. 

n=l 

In addition, in this case, we also have that 

^ oo 

inf Iy{^) = l3,{V)--Y,n^n{Vr. 

MSP -2,2 2 ^ 

n=l 

Proof. Condition (50) and the previous Proposition guarantee that there is a measure with 
full support such that J^^ Tn (|) fi{dx) = — 

The other way around works as follows. Assume that /ly is the equilibrium measure on 
[—2, 2] and has full support. What we need to show is that (50) is satisfied. 

Let an = /^2-^'i (f) f^vidx). Then, for any other measure u with /y(z^) < oo on [—2,2], 
from (45) we obtain that 

ly ((1 - e)^^v + ev) = Po{V) ^^r.{Vf + ^ ((1 - e)a„ + 6< + ' 

n=l n=l ^ ^ 

where a'^ = f Tn (|) u{dx). Since Jy(/iy) is the minimum of Iv{^) over all probability 
measures on [—2,2], differentiation with respect to e > at yields 

(51) Z^~("" + 2 J -0- 

n>l ^ ^ 

Now we consider measures of the form 

v^dx) = (1 + 5(f){x))^v{.dx) 
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where is a polynomial such that J (pdfi = and 6 is small in absolute value. Applying this 
for ±5, in (51), where S is small enough, we obtain that 

nf^nivy 



(52) 



^ n 

n>l 



T 



X 



(f){x)fividx) = 0. 



A word of caution is in order here. We need to justify that the measure u has finite loga- 
rithmic energy, namely that 

2 



n>l 



T 



X 



(1 + 5(t){x))iJiv{dx) ) < oo. 



This actually follows easily for each polynomial cj) = Tk for some /c > from the fact that 
2TkTi = T\k_i\ + Tk+i for any k,l>0. 

Because of (31) we have that 2 /log \x — y\fiv{dy) = V{x) + C, almost surely (with respect 
to the Lebesgue measure) on [—2, 2], and from Theorem 4.1, the density g{x) of fiy is given 
by 



g{x) = AiV4: - x^ 



V'{y) - V'{x) 



-dy + 



A2 + Asx 



-2 a/4 - _ a;) " a/4 - 

for some constants Ai and A2. In particular, since V is it implies that g{x) = ^^y^z^ for 
some continuous function h. Since the measure fiy has full support, h{x) > on a dense set. 

dx 



Next, we observe that because (3n{V) is square summable in {l[_2,2]{x)- 



R(x) 



E- 

^ n 



n>l 



an + 



nf3n{V) 



T 



is convergent in L^(l[_2,2] (a; 



dx 



therefore we deduce from (52) that 



R{x)(f>{x)h{x) 



dx 



V4 



0, 



x^ 



for any polynomial (p. This easily implies that R{x) = almost everywhere and this in turn 
results with a„ = —nf3n(y)/2. The rest follows. □ 

Theorem 4.2. // the equilibrium measure of a V which is and piecewise on [—2, 2] 
has full support then, 



'([-2,2]) '"^^ J^2nV4^ Jo \J^2 2nV4^J 



V{0)- 



'1 

t 



txV'{tx)dx 
27rV4 - x2 

Proof. According to Corollary 4.2, we have 



tV'{tx)dx 



2 7rA/4 — x"^ 
Iv = f5o{V)-\Y,nnvf. 



^ V'{tx)dx 
2 7r-\/4 — 

dt. 



dt 



n>l 
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Thus our task is to prove that 

Iv- n ,.r^2 r \ f r xV'{tx)dxV f V'{tx)dx\ , 

A polarization argument shows that this is equivalent to proving that for any potentials 
Vi and V2 

n ,.r^n /t.n /"^ ( C xV[{tx)dx\ ( xV;,{tx)dx\ , 



-2 n\/4: — x^ J \J-2 7rV4 — x^ , 

To do this, because of the linearity in Vi and V2 and the fact that polynomials are dense 
(with respect to topology) in the set of smooth functions on [—2,2], it suffices to check 
this for Vi{x) = x^ and V2{x) = x"^. If /c or m is zero, both sides of (53) are zero, therefore 
we need to check this for /c, m > 1. 
Now, we use 

= ^ E + ^ (';) . and = ± (*' + T,„«...(.) 

fc=0 ^ ^ ^ A;=0 ^ ^ 

from which a direct calculation yields 



(54) / x'T, 



2 



2/ 7rv/4^^ 



i—n I ' 
2 



with the convention that (p_,_\/2) ~ ^ p G Z and (*) = for p < 0. Therefore, (53) 
becomes in this case 



m \ / k \ mk / m\ / k\ mk fm — l\ fk 



2 ■ A^y V^y 4(m + k) Vf ; vi; (m + a;) v ;v 

In the case m, k have different parities, then the above expression is 0. If they have the same 
parities, the equality follows from the next Lemma. □ 

Lemma 4.3. The following identities hold 

2^1 \ / 2/2 \ ^1^2 / 2/i\ / 2/2 



J2p 



h - pj Kk - pj 2(/i + l2)\hj\ k 
2/1 + 1\ /2/2 + l\ (2/1 + l)(2/2 + 1) [2k\ [2I2 



h - P J Kh - P J Zi + ^2 + 1 \h J \ h 

V 

with the convention that (;^) = for g < or g > j. 

Proof. These identities can be checked with the zb package written for Mathematica. For 
details on this we refer the reader to the wonderful book [PWZ96]. For completeness we give 
here the main calculation. 
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The first identity is equivalent to 
Let us denote 

The idea of the zb method for our case is to write 

(55) fih + 1, l2,p) - f{lu h,p) = 9{lu h,P + 1) - 9{lu h.p)- 

and this proves that for all /i > 1 one has h{li,l2) = h{l,l2)- Since it is is immediate to 
show that h{l, l^) = 1, the rest follows as soon as we know that g{li, l2,p) = for p = 1 and 
for large p. 

The whole point of the zb method is to actually compute the function g{li, hiP)- We will 
refer the reader for the details to [PWZ96] and will give here just the results obtained with 
Mathematica. 

2p(p-i) (?.■;') CK) 



9{h,h,p) 



Notice that for p > min{/i + 2,l2 + 1}, g{hJ2,p) = 0. One can directly check (55), by 
dividing both sides by f{li,l2,p) which reduces it to an identity in the field Q(/i,/2,p)- 
For the second identity, as in the preceding argument, we want to show that 

^(2p+l)(/i + /2 + l)('^+^) (f+^) 



Defining 



f{h,k,p) 



(2/i + l)(2/2 + l)(^/;)ffl 

(2p+i)(/,+/2+i)eK)a-;) 



(2/i + l)(2/2 + l)R;^)(t^) 
the corresponding companion in this case is 



9{h,l2,p) 



(?,'«)(?) 



Equation (55) is satisfied and g{h,l2,0) = and g{h,l2,p) = for p > min{/i + 1,^2 + !}• 
This proves that h{li,l2) = h{0,l2)- Now, /i(0, ^2) = 1 which ends the proof. □ 

Before we state the next result, for a potential V, we define 

"2 V'{cx + b)- V'{cy + b) dy 



(56) i^bA^) ■-- 



x-y tiJA - y2 
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Theorem 4.3. Assume V is an admissible potential on R. Then the equilibrium measure on 
M associated to V has support the interval [—2c + 6, 2c + b] if and only if (c, b) is the unique 
absolute maximizer in Ml x M o/ 



(57) 



1 /"^ dx 
H{c,b) ■=\ogc- - V{cx + b)^= 



x^ 



and 

(58) Vb.c > on a dense subset of [—2,2]. 

If in addition V is a C"^ and piecewise potential on a neighborhood of the support [—2c + 
b, 2c + b], then [b, c) is a solution of 

' cxV (cx + b) = 2 



(59) 



J-2ncx + b)^^^=0. 



In this case the equilibrium measure fiy is given by 



and 
(60) 



— log c 



Hv{dx) = l[-2c+b,2c+b]{x) 



^ V{cx + b)dx 
2 7rV4 — x"^ 



^b,c{,{,x - b)/c)^Ac^ -{x-bf 
2c7r 



dx 



Vib)-\ogc 



-1+1 



^ xV'{sx + b)dx 

2 2TTy/4 - x^ 

sxV'{sx + b)dx 



27ly/4 



x^ 



+ 



^ V'{sx + b)dx 
2 7rV4 — x"^ 

^ sV'{sx + b)dx 

2 TT 



ds 



x^ 



ds. 



Proof. If the support of /iy is the interval [—2c + 6, 2c + 6], we have to prove first that (c, b) 
is the unique absolute maximizer of H . The function H appears in the literature as the 
F-functional of Mhaskar and Saff (see for instance [ST97, page 194]) and for the sake of 
completeness we adapt the proof of this first part from there. 
Define the arcsine law of the interval [—2c + 6, 2c + b] to be 

dx 

UJc,b{dx) = ll-2c+b,2c+b]{x)- 



v/4c2- (a:-6)2' 



A simple rescaling of equation (37), shows that / log \x — y\ujc,b{dy) > log(c) for all x, with 
equality only for x G [—2c + b,2c + b]. 

Integrating equation (31) against the measure Uc'^b' yields that 

log \x - y\uJc',b'{dx)fividy) > C 



2 logic') > C. 




V{x)ujc',b'{dx) — 2 
and thus, after interchanging the integrations 

V{x)uc',b'idx) 
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Because (31) is equality quasi-everywhere on [—2c + b,2c + b], this implies that we have 
equality in the above inequality for c' = c and b' = b. In fact, this is the only case of equality 
as otherwise 

C = j V{x)uc',b'idx) -2\og{c) > j V{x)uJc'^h'{dx) -2 j j \og\x - y\uc',b'idx)iJ,v{dy) > C, 

hence we must have that /ly almost surely, log(c') = / log \x — y\u!c',b'{,dy), which according 
to (37) is possible if and only if Uc'^b' is actually equal to cu^b, or c' = c and b' = b. 
From (57), upon differentiation with respect to c and b, we deduce that 

(61) / cxV'icx + b) =2 and / V'{cx + b)— ^ = 

J-2 7rv4 — J-2 7rv4 — 

which combined with (43) proves (58). 

To prove the converse, notice that because (c, fe) is a maximizer of H, we have (61). It 
is then clear that the nv solves equation (43). What we have to prove is that this measure 
satisfies condition (31). To this end, it is sufficient to prove that for any 6' G M and c' > 
one has 

V{x) - 2 y log |x - y\fiv{dy)^ Uc'^'idx) > C. 

Switching the integration in the double integral, and performing some elementary calcula- 
tions, this inequality becomes equivalent to 

V{dx + b')dx 



Tly/A — X^ 




log \x - y\uJc',b'{dx)iJ,v{dy) > C. 



This inequality is equality for c' = c and b' = b, and thus C = —H{c,b). If c' and b' are 
arbitrary, the inequality is a consequence of the fact that the left hand side of this inequality 
is greater than or equal to —H{c', b') which in turn is by the hypothesis > —H[c, b). 

Identity (60) follows from Theorem 4.2 applied to V{x) = V{cx + b). □ 

In the case of even potentials, we know that the equilibrium measure is symmetric and 
thus in the preceding result we can always assume that 6 = and this deserves a special 
statement because of its simplicity. 

Corollary 4.4. If y is a C^, piecewise and even satisfying (5), its equilibrium measure 
is supported on the interval [—2c, 2c] if and only if c > is the unique maximizer of 

'"^ V{cx)dx 



H{c) = log c — / 
Jo 



7r-\/4 — 
and 

~ V'{cx) - V'{cy) dy 
-2 x-y 71^/4 - ?/2 
is positive on a dense set of [—2, 2]. In particular c solves 
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In this case the planar hmit is 

dx 



(63) 



Iv = V{0)-\ogc- T- -1+fl-/ sxV'{sx)^L==] 
Jo s \ Jo 7rV4-xV 



ds. 



We point out here an interesting property, namely, that the solutions (c, b) of the system 
(59) are critical points of the functional ly from (60). 

Proposition 4.5. Let V be a potential on M and consider 



Iv{u, v) = —logu + 
If{c,b) satisfy 



V{ux + v)dx 

■K\/4 — 



xV'{sx + v)dx 
27rV4 - x2 



+ 



2 T// 



y (sx + f )dx 



(64) 

If {c,b) satisfy (59) , i/ien 
(65) 



2 r(cx + &)t/x 
2 7ri/4 — x"^ 

div 



dv 



dh 



0. 



(c,6) 



(9u 



0. 



(c,fe) 



/n particular the critical points of H from (57) are also critical points of ly- 
Proof. To see (64), after differentiating with respect to v, we need to show that 



dl 



V 



dv 



icb) 



2 T// 



V'{cx + b)dx 
2 7r-\/4 — 



c /.2 /.2 



s(xy + 4:)V'{sx + 6)V^"(sy + 6) 
2 J -2 47rV(4-a:2)(4-?/2) 



dxdyds = 0. 



Now we present the following result. 



Lemma 4.6. If f/ G C^([— 2, 2]) or is a formal power series, then the following holds 
(66) 



2 47rV(4-a;2)(4-?/2) 



In particular, if [/ satisfies. 



f/(x) 



2 TT 



2 TTa/T 



= 0, 



x^ 



2 7rA/4 - y^ 



x^ 



then 



1 /•2 r2 



2 J -2 



s{xy + 4)f/(sx)f/'(s2/) 
47r2^(4-x2)(4-?/2) 



dxdyds = 0. 
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Proof. By polarization, it suffices to show that for any two potentials or formal power 
series, Ui and U2 on [—2,2], we have that 



sixy + A)iUi{sx)U^{sy) + U2isx)U[{sy) 



2 47rV(4-a;2)(4-?/2 



2 TT 



zdx 



x^ 



2 TV J A- y^ 



dy ds + 



dx dy ds 



U2(X) 



2 vrV4 



-.dx 



x^ 



sU[{sy) 



■.dy ds. 



It is clear now that it suffices to check this for Ui{x) = x^ and U2{x) = x 
help of (54), becomes 



2 71^/4:- y^ 

which, with the 



n + m 



n n 



4 \ - 

^ ^ 2 



m 



m + 1 



m + J 

■m+l 
2 



m I m 



n 



71+1 

2 



+ m 



m — 1 

m— 1 



n 



n + 1 



n 

n-l 
2 



m — 1 

m— 1 



+ n 



n — J 

11-1 
2 



Here we use the convention that (^) = if 6 is not a nonnegative integer. As long as n and 
m have the same parity, both sides of the above expression are 0. Also due to the symmetry 
in n and m, it suffices to check this for n = 2k and m = 2/ + 1. In this case, it is easy to 
prove that both sides are equal to 



21 + 1 f2k 



2l + 2\k 



□ 



Taking U(x) = V'{cx + h) in the Lemma, after a simple change of variables, the rest 
follows. 

Equation (65) is clear from the fact that 



/y (u, v) = V{v) — log u 



"1 
s 



-1+1 



^ sxV'{sx + v)dx 
-2 27r\/4 - x2 



+ 



^ sV'{sx + v)dx 
-2 7r\/ 4 — x-^ 



and thus the u-derivative vanishes under the condition of (59). 



ds. 
□ 



5. Examples and Computations with Analytic Matrix Models 

5.1. Cases of One-Cut Potentials. With the result from Theorem 4.3, it is instructive 
to recover the classical results (see [ST97] where weaker regularity conditions are required) 
which guarantee that there is a one interval support of the equilibrium measure. 

Corollary 5.1. Assume that a potential V satisfying (5) is either convex or even with 
xV'[x) increasing on [0, 00). Then the equilibrium measure has one interval support and 
the maximizer is non-degenerate (i.e. the Hessian of H is invertible at the maximizer). In 
addition, the function iIjc^, is positive on [—2,2]. 

Proof. First, we need to check that the function H{b, c) has a unique maxima. 
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In the case V is convex, we show that H is concave. Indeed, the hessian of H at (c, b) is 



(Hessi7)(c, b) 



and strict concavity is equivalent to 



^ x^V"{cx + b)dx 

2 27rV4 - a;2 
(ex + b)dx 

27rv/4 - x2 



2 27rV4^^ 
2 27rv/4^^ 



1 



> and 



V"{cx + 6) da; 



2 27rV4 - 

^ V"{cx + /"^ xV'(cx + b)dx 



27rV4" 



x^ 



27rV4" 



x^ 



27rV4 - x2 



27rV4 - a;2 



> 0. 



The only way either of these fail is if V"{cx + 6) = for all x G [—2,2], which implies 
V{x) = Ax + B for some constants A, B and all x G [—2c + 6, 2c + 6]. This in turn results 
with F{c',b') = logo' — B for all c' < c which contradicts the assumption that (c, 6) is a 
maximizer of F. 

On the other hand, one can easily check that H is concave on (0, oo) x M. This combined 
with strict concavity near (c, b) implies that the maximizer is unique. 

In the case V is even and xV'{x) is increasing, we may assume that 6 = and thus the 
function H becomes a function of one variable with 

xV'{cx)dx if cxV'{cx)dx^ 
c 



H'(c) 



1 

c 



27rV4 



x^ 



vr 



x^ 



Now, since the function xV'{x), is increasing, one can see that cH'{c) is decreasing and thus 
there is only at most one critical point of H. On the other hand, one can check that there 
is a maximizer of H{c), hence we deduce that there is a unique such maximizer. 
In addition to this, the Hessian of H{c, b) at the maximizer (c, 0) is 



(Hessi/)(c, b) 



2 x^V"{cx)dx 
-2 27rA/4 - a;2 











2 V"{cx)dx 
2 27rv^ 



x^ 



Now, using the fact that H'{c) = and a simple integration by parts reveals that 

"2 c^V"{cx)V^^dx A^V"{cx)dx 

-2 



1 



^ cxV'{cx)dx 
.2 27ryr 



x^ 



27r 



2 27rV4^ 



2 c^x^V"{cx)dx 
-2 27rV4 - x2 



which implies 

(*) 



V"{cx)dx 



1 



_2 27rV4 - x2 c2 



2 a;2l^"(cx)(ix 
2 27rV4 - x2 ' 



Now, if we denote g{x) = xV'{x), then g is an increasing function on [0, oo] and therefore 
x'^V"{x) = xg'{x) — g{x) > —g'{x) for all x > 0. In particular we obtain that c2x2V^"(cx) > 
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-g{cx) for all x G [0,2]. Furthermore, from the equation determining c, we get 

g{cx)dx 



1 



7r-\/4 — 



which in turn implies 

1 , c^x'^V"{cx)dx ^1 1 /"^ ^(ca;)rfa; _ ^ 



Jo 7rV4 - x2 c2 c2 7o 7rA/4 - 
and this means that quantities in (*) are positive, thus the Hessian of F at (c, 0) is non- 
degenerate. 

Having checked the uniqueness of the maximizer, we need to check the other condi- 
tion. In the case of convex potentials, the non-negativity of ipcb follows from the fact that 
V {cx+b)-v [cy+b) y g x,y E [—2,2]. Furthermore, t/Jcfiix) = 0, enforces V'{cx + b) = 

V'{cy + h) for all y G [—2, 2], which in turn yields V{c ■ +h) is constant on [—2, 2], something 
which is contradicted by the assumption that (c, 6) is a maximizer of H{c,b). Hence we 
actually obtain the stronger conclusion, namely ipc,b{x) > on [—2,2]. 

In the case V is even and xV'{x) increases on [0,oo], one can show that ipc is an even 
function and with simple manipulations of integrals that 

'"^ xV'{cx) — yV'{cy) dy 



ijjcix) 



which makes clear that ipc{x) > for all x G [—2, 2]. □ 

5.2. Analytic planar limits of various even potentials. In this section we explicitly 
compute the planar limit of some 1-cut potentials, illustrating the formulas of Section 4. A 
typical example is the case where V is a smooth potential which is analytic near the support 
of the equilibrium measure. 

The easiest to deal with is the case of even potentials because in this case we can invoke 
Corollary 4.4 and reduce the problem of determining the support of the equilibrium measure 
to the maximization of a function of a single variable. In this case the planar limit is actually 
a one variable function of the right endpoint 2c of the equilibrium measure. 

Assume that V is an even potential such that it has a power series expansion valid on a 
neighborhood of the support: 

oo 2n 

(67) V{x) = J2a2n^. 

n=l 

In this case, from Corollary 4.4 we learn that 

, , 1 >^ a2nc2" x^^dx , 1 ^ a2„c2" f2n\ 
ii c = log c > / = log c > 

n=l " ^ V n=l ^ ' 

where in the last equality we used equation (54). The critical points of this function satisfy 
(62) which becomes 



(68) E«-c^"(t)=2- 

n=l 
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If c is the maximizer of F, then, again from Corollary 4.4 and (54), the planar limit is given 
by 

2' 

-1 



- log c + 



(.4£...-r:)) 



dt. 



Example 5.2. For V{x) = a2n^, with a2n > 0, and n > 1, the support of the equilibrium 
measure is [—2c, 2c], where 



1 

' 2n 



1=0 



a2n 1 2n 
. 2 U 

In this case, the equilibrium measure is 

/iy(f/x) = l[-2c,2c](x)— ^c(x/c) V4c2 - x^dx, ^,(x) = aa^c^-^ ^ ( jx^^^-'-D 

;— n \ / 

and the planar limit is 

j^_loga2n , log(r:)/2) ^ 3_ 
2n 2n 4n 

To see this, one has to realize that (68) becomes in this case 



n ^ 

which has only one positive solution, and this is the maximizer of H(c) = logc — °^2w^" Cn) • 
The rest of the equalities are straightforward calculations. 

It is worth pointing out that in this example the potential is convex and thus, the equi- 
librium measure must be supported on a single interval. 
For n = 1, we recover the semicircular law. 

Example 5.3. Assume V{x) = a2n^ + ^2™^^ with > and 1 < n < m. In this case 
the equilibrium measure has a single interval support if and only if 

(69) a2n > —CnmO'2rn 

where 

m — n 

In this case, the support of fiy is [—2c, 2c] where c is the unique positive solution to 
(71) a2nc'-(^'')+a2mc'-(^"')=2, 
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the equilibrium measure is 

ljiy[dx) = ^^l[^2c,2c]{x)iic{x / c)y/Ad^^^dx 

n— 1 /r\i\ 1 

,2(m-/-l) 



1=0 ^ ^ 1=0 ^ ^ 

and the planar limit is 



X 



logc+-^ + 



2n V n / 2m \ m 



4(n + m) \n J ym J 16n \ n y 16m \ m 
To prove these, we need to look at the critical equation (68) and notice that for 

/(C) = a^nC^" V ^^n.C^'" V 2, 



one has 

It is clear that /' has at most one positive root. If cq > is the positive root of /', then 
/'(c) < for < c < Co and /'(c) > for c > cq. If /' does not have any positive root, then 
/' > 0. Since /(O) = —2 and /(oo) = oo, it follows that / must have a unique zero which in 
turn is the unique maxima of H{c). 

Now having proved that there is a unique maxima, we need to check the second condition 
from (4.4). That boils down to 

Mx) > 

on [—2, 2] with strict inequality on a dense set. This is equivalent to 

v^m-l (2l\ 2{m-l-l) 
^ _ 2m-2n Z^i=0_U/2 

Er=o ( /)^'^"-'-'^ 

for all X G [—2, 2] which in turn is satisfied if and only if 

where Knm is defined by (70). On the other hand from the critical equation (71) replacing 
a2n, we arrive at 

2-a2mc2'"(''") 

/2n\ — '^2mC J^nm 

\ n ) 

and thus, after noting that Knm < (^™) / (^^) , is the same as 

l/{2m) 

/ 2 \ 



a2m{{Z) - ^nm(^)) 
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For the function /, we know that f{x) < if and only if a; < c. Thus, we have the second 
condition in Corollary 4.4 satisfied if and only if 

which is equivalent to equation (69). 

The constant Knm from (70) depends only on n and m. It can be explicitly computed in 
the case n = 1 and any m > 2 as the minimizer is t = and thus Kim = {^^-2) I ( nT-f ) 
then a simple rearrangement reveals that 

l2m-2\ 
\ m-l ) 



Clr 



/2m-2\ ' 
\ m J 

In general, it does not seem that one can find an explicit algebraic expression of the minimizer 
in (70). For the case of n = 2 and m = 3, we have an exact solution as the minimizer in the 
expression there is t = —2 + \/6 and then in this case K23 = 2\/6 — 2 which produces 



C23 = V 4 + v^- 

The root c from equation (71) does not have a simple representation in general. However, 
in some cases it can be solved explicitly. For example if m = 2n, one has 



and similarly there are algebraic expressions in the case m = 3n or m = 3k, n = 2k and also 
m = 4n or m = Ak, n = 3k, but we omit the lengthy formulae here. 

Corollary 5.4. For the quartic potential 

V{x) = flay + 

the equilibrium measure has a single interval support if and only if a2 > —2^Jai in which 
case 



h{x) = ^l[^2c,2c]{x)V4:C^ - x^{b2 + a4X^) 



-0.2 + V '^2 + 12a4 
6a4 



b2 
W 



2a2 + \/al + 12a4 




860204 + 162a| + {al + 800204) y^a^ + 12a4 



2 / 482a| 



We should point out that this example appears for instance in [Joh98] 
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6. Matching Formal and Analytic Matrix Models 

In this section we will prove Theorem 1.2. Our first task is to match the analytic equations 
Equations (59) of (6, c) of a 1-cut potential with the equations (13) for {71, S). Consider a 
1-cut potential V and its Taylor series expansion at x = 0: 



V{x) = 



n 

n=l 



Using the key identity 
(72) 



x'^dx j (^/g) if is even 



-2 7rA/4 — 1 if n is odd 

and interchanging summation and integration, (54) gives 



n>l j>l 



Then, Equation (59) gives the system of non-linear equations for [b, c) 



(73) 



Following the notation of [BDFG02], let us use the change of variables (6, c^) = {S, R) as in 
Equation (8). Then, (i?, S) satisfy the system of equations 



(74) 



Consider now the 1-cut potential 



n=l 

Then, Equation (74) gives the system of non-linear equations for {R, S) 



(75) 

Using 



it follows that {R, S) satisfy the system of non-linear equations (13). 

Observe that for a fixed admissible potential V, Equation (13) may have none or more 
than one real solutions for {R, S) but for small parameters a = (ai, 02, ... ) in some £,J for 
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small enough r, R and S become analytic functions of a (see Theorem 12.1). However, it 
always has a unique formal solution (7l,S) G (1 + A~^,A^). 
This proves that TZ = R and S = S in Theorem 1.2. 

To finish the proof of Theorem 1.2, we need to prove that the coefficient of any monomial 
a"^ . . . a^'^ from the power series J-q and Fq are equal. The important point here is the fact 
that the each such monomial involves finitely many ai, 02, . . . , and thus we may assume 
that all the 1-cut potential is actually a polynomial. 

Now, assume that all for n > k and consider potentials of the form 



\n=l 

For small real parameters ai,a2, ■ ■ ■ the functions 



n 2k + 2 



1 , /^^exp(-iVTr(V(M)))rfM 



g^{a,, a„...,a,)-— log ,^p(_;vTr(MV2)))dM 

are analytic in ai, . . . , on a neighborhood of G M'^. Since the limit g^o exists, the limit 
is going to be also an analytic function in these variables. This means that at the level of 
power series the coefficients must converge to the coefficients of the limit. 

On one hand expanding the (^at in power series, the limiting coefficient of a"^ . . .a'^'' is 
exactly the corresponding coefficient from the formal model. On the other hand, the limiting 
function g^o is obtained via the potential theory and using the perturbation theory from 
Section 12, it is easy to see that c, b, the solution of the system (75) are actually analytic 
functions of ai, . . . , a^. In particular it means that the planar limit Fq is equal to g^o from 
(28) and is analytic, thus concluding the proof. □ 

Remark 6.1. From now on, whenever we have a formal potential V = ^ — o,n^, we 

will use b, c as the solution to (73) which has a unique solution in (c, 6) G (1 + A"^, A~^). 

7. The planar limit J'o(^) in terms of TZ{t) and S{t) 

This section is devoted to the proofs of Theorems 1.3 and 1.4. After we discuss the 
proofs we give a main consequences of these formulae, namely the fact that the planar limit 
enjoys algebricity in some cases which allows complete description of the asymptotics of the 
coefficients of J-q. 

7.1. Proof of Theorem 1.3. In this Section we will prove Theorem 1.3 and the first part 
of Remark 1.3. In this section, it will be convenient to use the 1-cut potentials Ve and Ve 
given by 

2 r) 

/'I ry''' 



n>l 



^ n 

n>l 
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For simplicity of notation in this section we will drop the dependence on e from the writing 
of Ve and Vg. 

We start by setting c{t),b{t), and c{t),b{t) to be the power series solutions to (73) corre- 
sponding to potentials V, respectively V. From the fact that V(t,x) = V{t, \/ix), we easily 
get that 



(76) 

Then V(t, x) 

(77) 



X 

2t 



c{t) = yftcit) and h{t) = Vtb{t). 
— yV{x) and the system satisfied by c(t),b(t) is given by 

<2 -/ 



J^^~c{t)xV{trc{t)x + b{t))^^ 
J^^V{t,c{t)x + b{t))-^ = 0. 



where the derivative V'{t,x) is taken with respect to x. Set now, 



X(t) = -log£(t) 



V{t,c{t)x + b{t))dx 



2 

m 



IT 



X^ 







xV{t,sx + b{t))dx 
27rV4 - x2 



2 -r,, 



V'(t, sx + b{t))dx 



IT 



ds. 



Taking the derivative with respect to t, 

c'{t) {d'{t)x + b{t)))V{t, d{t)x + b{t))dx /"^ V{t, 5{t)x + b{t))dx 



m 

- c'mt) 

- 2b' {t) 



X^ 



IT 



X^ 



xV{t,c{t)x + b{t))dx 



x^ 



2 -Tv 



V{t,c{t)x + b{t))dx 



TT 



c(t) 1-2 1-2 







s{xy + 4)V'(t, sx + b{t))V"{t, sy + b{t)) 



^2J~2 47rV(4-x2)(4-?/2) 
/-s /-^ s(xy + 4)V'(t,sx + 6(t))l)'(t,s2/ + 6(t)) 

2 J-2 



dy ds 



47rV(4-a;2)(4-?/2) 



(is, 



where V(t,a;) is the derivative with respect to t. Since V(t, x) = — the system (77) and 
Lemma 4.6, we can simplify this to 



X'it) 



1 
2 



2 /- 



c(t)x + b{t)ydx 



TT 



c(t) /.2 /.2 



X^ 



s(x?/ + 4)V'(t, sx + b{t)){sy + 6(t)) 



J-2J-2 47rV(4-a;2)(4-?/2) 
25(t)2 + 6(t)2 ^ 1 r^*^ r g(g^ + 2&(t))V^(t, sx + b(t)) 
^^Jo J -2 7rv/4^ 



2^2 



(ix (is 
(ix rfs. 
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Next, observe that for any continuous function / : [—2c, 2c] — )■ M with c > 0, one has 



Fubini J-2Jc\y\/2 IT W 4s'^ - C^y"^ 



^2 (ci/ + 26)/(cy)v/4^ ^ 
7-2 47r ^ 

= -/ (z + b)fiz-b)V^c^-iz-b)Mz. 

y={z-b)/c 47r y_2c+b 



Going back to the previous equation we now have 
•■m s{sx + 2b{t))V{t, sx + b{t)) 







-2 7r-v/4 — 

471" J~2c{t)+b(t) 

Take the derivative with respect to t and observe 

^ l'2c{t)+b{t) 
-2c{t)+b{t) 



dx ds 

+ kt))V{t, z) \]Ac{tf - {z-h{t)Ydz. 



dt 



[z + h{t))V{t, z)yjAc{tf - {z-b{t)fdz 



2c{t)+b{t) 



2c{t)+b{t) 



b'{t)V{t, z)^Adity - {z-b{t)fdz 



^2 



2c(t)+fe(i) ' y/4£(t)2 - {z-b{t)y 

{z + 6(t))"i>'(t, 2) Y^4g(t)2 - (z-6(t))2ci^ 



2c{t)+fe(t) 



2c(t)+6{t) 



2c(t)+6(f) -2b{tfh'{t) + 4c(t)26'(t) + 4g'(t)c(t)6(t) + z 2b'{t)b{t) + 4£'(t)c(t) 

V'(t, 2) ^ '-dz 

-2cit)+m \kc{tY -{z- h{t)y 



2c{t)+b{t) 



{z + my Ac{tY - {z -b{t)Ydz. 



2c{t)+b{t) t 

Changing the variable z = c{t)x + h{t) and using the system (77), we obtain 

J l'2c{t)+b{t) ^ _ I 

{z + h{t))V'{t,z)JAc{tY-{z-b{t)ydz = 

-2c{t)+b{t) 

A7r{b{t)b'{t) + 2d{t)d'{t)) - 2nd{tf{2bitf + 



dt 
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Therefore we arrive at the equation 



{ex\t))' = - -c{tf - ^ + + 2c(t)c'(t)) - 



dt\ ^ ' 2 / ' ' ' ' ' y > y 

C(t)2(26(t)2 + g(t)2) 



2^2 

Since J-o,e(^) = I —^(t), it imphes 



(t^-^0,e(t))' 



27^e(^)52(^) + 7^2(^) 



2 

This is exactly the statement from (20). To prove also the statement from (19), namely that 

1 /•* (t-s)(27^e(s)5,^(s)+7^g(s)-l) 

denote the right hand side by Q{t) and notice that both sides satisfy the same differential 
equation, namely 

In addition, a direct check reveals that 

Ae(0) = 6^(0) = 0, J-^,,(0) = 6;'(0) = a2/2 + a2/2 

which actually follows from the fact that 7?.(t) = 1 + ait + 0{t^) and S{t) = y/tai + 0{t) 
(see for example the formulae in Appendix A). 

7.2. Proof of Theorem 1.4. In this Section we will prove Theorem 1.4 and the last part 
of Remark 1.3. It will be convenient to use the 1-cut potentials V/ and V/ given by 



2 ^ n 



n>3 

k 



k>3 

As we did in the previous section, for the sake of simplicity we will drop the dependence 
on / from the notation V/ and V/. 

Define c{t),b{t) and c{t),b{t) the power series solutions to (73) corresponding to V and 
V/t. Then, one can easily check that 

(78) c{t) = Vtc{t), b{t) = Vtb{t). 

The corresponding system of equations for c(t) and b{t) is 



(79) 



r_,~cit)xVicit)x + bit))^^ = 2t 
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Now set ^o(^) = "'^^^ov/f Thus 



Goit) = t'\ogc{t)-t 



V{c{t)x + h{t))dx 



TX 



xV{sx + b{t))dx 



x^ 



"2 -V,/ 



V{sx + h{t))dx 



L2 tta/I^ 



x^ 



ds. 



Differentiating this with respect to t and keeping in mind the system (79), we get 
Goit) = 2t log c{t) - 



2 7rV4 — x"^ 

+26'(«) r r f ^(^■■i/+4)v-(..+Mt))K(.,+sw),. 

Jo J-2J-2 4„^-y'(4-x^)(4-y^) 

Now, taking U{x) = V'{c(t)x + b{t)) in Lemma 4.6 and a simple change of variables proves 
that the last term of (80) becomes 0. Thus, we can continue (80) with 



V{c{t)x + h{t))dx 



2 7rA/4 — 



Differentiating this with respect to t, and using again the equations from (79), we obtain 

G'^{t) = 2\og c{t)- 

In other words, integrating this twice and keeping in mind that ^o(O) = Q'q{^) = 0, we get 



Goit) = 2 [ {t-u)\ogc{u)du. 
Jo 



Now, one has to notice that an easy calculation yields, 

3^2 ^2 rt nt 

t^Toit) = — - -\ogt + Go{t) = 2 / {t - u) log c{u)du = / log7^/(n)c^u, 
4 ^ Jo Jo 



from which Theorem 1.4 follows. 



7.3. Proof of Proposition 1.8. Part (a) and (b) of Proposition 1.8 follows from Theorems 
1.3 and 1.3. 

Part (c) and (d) follow from Proposition 1.7. 



8. The planar limit for extreme edge potentials 
In this section we will compute the planar limit for five extreme formal potentials. 
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8.1. Exact Formulae. Consider the extreme formal potentials Ve{x), V°^{x) E Q[[[t^^'^]][[x]] 
given by 

(80) Vr(x) = ^ + ^log(l-tx^) = ^-$:*^ 



2 2°' ' 2 ^ 2n 

n=l 

(81) Ve(x) = y+log(l-v^x) = ^-5^^ 

n>l 

These potentials correspond to counting planar diagrams with even respectively arbitrary 
valency of the vertices and a fixed number of edges. Their corresponding invariants h = 
= b{t) = Se{t) is an element of Q[[t^/^]], while c = c{t) while Se = Se{t), 7^e = 7^e(^) 
and J-Q^e = J^o,e(t) are elements of Q[[t]]. Our next proposition summarizes the algebraic 
properties of these elements. 

Remark 8.1. For simplicity of writing, in this section we will omit the subscript e in writing 

Proposition 8.2. (1) For the potential V|^, we have 
S{t) = 

1 + At- v/T^8t 



f82) 



Mt) = T7^2 o log 



t 

24t + 72^2 - (1 -20t)yr^8t 3,^^l-4t + yr^8t 
128t2 



1 + :dL + + + + 132^6 + 1^ + 3432^8 ^ ^(^lO) 

2 4 5 7 



(2) For i/ie potential Ve, we have 



(83) 



5(t) 

m 



m - VI - i2t 

18t 



^ , , 1 - 36t + 162^2 - (1 - 30t)Vl - 12t 1 , l-6t + ^1 - Ut 
^'^'^ = 216^^ 2^°^ 2 

189^4 1458^5 8019^6 104247t^ 

= t + — + 9t^ + ^— + + — — + + 0{t^) 

4 4 5 4 7 ^ ^ 

Proof. Solving the nonlinear system of Equations (13) for our formal potentials Ve(x) and 
Vg^(x) seems at first an impossible task. Instead, we will use the analytic ideas from Section 
4 to translate this system into a more tractable one. 

In Section 6 it was shown that Equations (13) for {b,c^) = {S,TZ) are exactly Equations 
(59) for [b, c) in case of admissible analytic potentials. The proof also works for formal po- 
tentials, too, such as our potentials Ve(x) and V°^(x). Proposition 4.5 implies that Equations 
(59) are the critical point equations for the function T-Lib^c) from Equation (57). The last 
function can be computed explicitly for the two formal potentials Ve{x) and V*'^(a;). 
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To prove part (1) of the proposition, is even so b{t) = 0. Computing the function "H 
gives 

1 /"^ log(l - tc^x^)dx 



If 

n{c) = logc- — - - / 



2 7r-\/4 — 
1 , 1 + VI - 4te2 



^2 



log C lo; 



2 2° 2 
and thus 

•H'(c) = - 2c2 + ^ 



2c V VI -4te2, 

The solution c to 'H'(c) = such that c(0) = 1 satisfies a quartic equation 

4c^t + c2(-l - 4t) + t + 1 = 

and it is given by 



c(t) 



1 + 4t - vT^St 



8t 

Given and c(t) together with (19), we get (82). 

For part (2) of the proposition, observe that By Equation (37) we have: 



,2 



n{b,c) = logc- — - — -- 



6^ 1 log(l - Vtcx - Vib)dx 



2 4 2 y_2 ttVT 



2 ^2 1 Vt5)2-4tc2 

log(c) log . 

6V ; 2 4 2 2 



The critical point [b, c) satisfies the system 



1-c' 



,2 I 2tc^ 



M) 



-6+ , ^ = 0. 



Solve for \/ (1 — Vtb)^ — Atc^ = y/i/b and plug it in the first equation which becomes 



1 - c^ + '—^ p — = 0. 

{Vi/b){l-bVi + Vt/b) 

In turn, this implies 

2 _ b + Vt-b'^Vi 

^ ~ b + Vt-Sb'^Vi' 
We need to pick the solution c which for t = is 1 and thus 

_ Vb + Vt-b^Vi 

Vb + Vt-3b^Vt' 

We go back to the second equation of (84) and solve for c as a function of b to get 



2cH 
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Equating now the two expressions of in terms of b shows that 

b + Vt-b^Vi _ (1 - bVifb^ - t 
b + Vi-3b^Vi~ 4:tb^ 

This imphes that b satisfies 

{-b-Vt + b^Vtf{-b + Vt + sb^Vt) = 0. 

There are four solutions to this equation, 

1 - VI - 12t 1 + VI - 12t l-VrT4t i + v/rT4t 

Since 6(0) = 0, this ehminates the second and the fourth solutions. To decide which one is 
the right one, we notice that c(0) = 1 and this implies 



1 - VI - 12t , /I + 12t - VI - 12t 
= ^ and c = \ . 

Now using (19) one concludes (83). □ 

8.2. A review of holonomic functions and their asymptotics. In this section we briefly 
review soem standard facts about holonomic functions and their asymptotics from [PWZ96]. 
Recall that a formal power series 



(85) f{x) = J2<^n^ 

n=0 

is holonomic if it satisfies a linear differential equation 

d 

5^c,(x)/(^)(x) = 

j=0 

where Cj{x) G Q[x] for j = 0, . . . , (i with q(x) 7^ 0. A sequence (a„) is holonomic if it satisfies 
a linear recursion 

r 

^lj{n)an+j = 

j=0 

for all G N where 7j(?t-) G Q,[n] with 7r(?T-) 7^ 0. It is easy to see that a sequence (a„) is 
holonomic if and only if the generating series (85) is holonomic. Of importance to us are 
algebraic functions y = f{x) i.e., solutions to polynomial equations 

d 

(86) J]c,(a;)y^ = 

i=o 

where Cj{x) G Q[x] for j = 0, . . . , c? and Cd{x) 7^ 0. Algebraic functions regular at x = are 
always holonomic; see for example [CC86]. Moreover, algebraic functions regular at x = 
have holomorphic extensions to a finite branched cover of the complex plane branched along 
a finite set of algebraic points, given by the roots of the discriminant of the polynomial (86) 
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with respect to y. Locally, at a point x = xq G Q of the field of algebraic numbers, an 
algebraic function y{x) has a convergent power series expansion of the form 



y{x) 



n=0 



for some natural number d and for algebraic numbers Cn/d- This is the content of Puiseux's 
theorem [BPR06, Wal78]. If y{x) is regular at x = with Taylor series 

oo 

y{x) = ^ctna;" 

n=0 

the asymptotics of the sequence (a^) can be computed explicitly by the singularities of 
y{x) which are nearest to x = 0. The computation also includes the Stokes constants. A 
computer implementation of the rigorous computation is available from [Kau09]. In fact, 
(an) is a sequence of Nilsson type discussed in detail in [Gar 11]. 

8.3. Holonomicity and asymptotics. In this section we illustrate Proposition 1.8 with 
the concrete examples of the extreme potentials and study the coefficients of the Taylor series 
[fn) of the planar limit written as 

oo 
n=l 

Proposition 8.3. (1) For the potential V|^, TZ, Qq = Go{t) = J^oit) and satisfy 

mH - 7^(l + 4t) + t + 1 = 
64^3^2 ^ ^4g^2 _ 24t + 2)^0 + 9t - 1 = 

3 + 6(4t - 1)J^Q + 2t(8t - 1) J'o = 0, J'o(O) = 0, J'o(O) = 1/2 
{n + 3){n + - 4n(l + 2n)/„ = 0, n > 1, A = 1/2. 

In addition, 

(m F(t)-\^ 3(2n-l)!2-^ 

n>l ^ ^ 

and for large n, the asymptotics of is 

3 8"/ 25 945 16275 / 1 
^ ' ~ 4v^n7/2 " 8^ + 128^ ~ 1024^ ^ 

(2) For the potential Ve, we have that S, 71, Qq = J-'q and J^q satisfy 

3ViS^ -S + Vt = 
9tn^ - (12t + l)7^+ 1 = 
(90) lOSt^g^ + (108^2 - 36t + 2)^0 + 27t - 2 = 

3 + 3(6t - 1) J'o + t {12t - 1) J-^' = 0, J'o(O) = 1, J^o(O) = 1 
(n + 3)in + l)/„+i - 6n{l + 2n)U = 0, A = 1. 
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In addition, 



, , 2(2n- 1)!3" „ 



n> 

and for large n, 

2 12" / 25 945 16275 / 1 
^ ^ ~ V 16^ ^ 256^ ~ 2048n3 + 

Proof. It is straightforward to see that (87) follows from (82) while (90) from (83). 

For (1), a direct check proves that -Fo(^) solves the third equation of (87). This immediately 
implies the recurrence on /„ and then the closed formula in (88) which in turn combined 
with Stirling's formula leads to (91). 

Another way of checking the closed formula (88) is the following. Observe by a direct 
calculation that 

(8^2 - t) J^"(t) + (28t - 4) j;^'(t) + I2r^{t) = 0, J-o(O) = 0, J-^(O) = 1/2, J^'(O) = 3/2. 
This is a hypergeometric equation, and its solution is given by 

, , I , , s 3(2n - l)!2"-i „ 

J-o(t) = -t 3^2(1, 1, 3/2; 2, 4; 8t) = \un'., C 

n>l ^' 

where 3-^2(01, 02, 03; 61, 62! x) stands for the hypergeometric function with parameters Oi, 02, 
and &i, 62- This is exactly (88). Using Stirling's formula (see [01v97]), one can easily deduce 
(88). 

For (2), use the same proof as for (1). □ 

8.4. Three more flavors of the extreme edge potentials. In this section we will inves- 
tigate the following three flavors of the extreme edge potentials (81) and (80), given by 

(93) V,(x) = ^^ + ^log(l-tx^) = y-E^ 

n>2 

(94) V2{x) = — + ^tx + \og{l-Vtx) = — -Y, 

n>2 ^ 

(95) V3(a;) = ^ — + v^a; + log(l - Vtx) = — - • 

^ ^ Til 

n>3 



These correspond to the counting of planar diagrams with a fixed number of edges and 
vertices of even valency greater or equal to 4, or arbitrary valency greater or equal to 2, 
respectively arbitrary valency greater or equal to 3. 
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Proposition 8.4. (1) For the potential Vi, we have 



(96) 


S{t) 


(97) 


m 


(98) 




(99) 




(100) 





1 + 5t- v^(l + t)(l -7t) 
8t(l +t) 

1 - 22t + - (1 - 19t) v/(l + t)(l - 7t) 1 
128^^ gMl + t) 



3 , 1 - 3t + J{l + t){l-7t) 

— ioff 

8 ^ 2 

5t^ 23t^ 51t^ 124^6 2515t^ 13245t^ „, 

= — + — + + + + + + 0(f). 

2 6 8 5 3 14 16 ^ ^ ^ 

(2) For the potential V2, we have 



(101) Sit) 

(102) 7^(^) 

(103) J-o(t) 



1 - 5t - Vl - 10t + t2 
6v^ 

1 + 14t + _ (1 + t)Vl - 10t + t2 
18t 

1 - 32t + 96t^ + 76^3 + t4 _ (1 _ _ 27^^ + t^)^! _ lot + ^2 

216^2 



, 1 + t + VI - lot + t2 
log 



/...N t 8^3 312^5 1076t6 15528t^ ....os ^/ 

(104) = - + — + — + m^ + —^ + ^- + ^^ + U508t' + Oit'). 



(3) For the potential V3, we have 
(105) S{t) 



I- At- Vl - 8t - 8^2 



(106) 7^(^) 

(107) J-o(t) 



6(l + t)v^ 
1 + 16t + 16^2 - (1 + 16t)Vl -8t-8t2 
18t(l + t)2 

1 - 28t + 6^2 + 176^3 + 142t^ - (1 - 24t - 78t2 - 52^3)^1-8^-8^2 

216t2(l +t)2 



/...N I1 N 1 1 + 2t + VI - 8t - 8^2 

(108) +- log(l + t) - log 

3^3 47^4 ^39^5 43Q^6 11175^7 75149^8 

(109) = — + — + + + + + + 0{t^). 

^ ' 2 2 8 5 314 16 ^^ 

Proof. We follow the same approach as in Proposition 8.2. 
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For (1), the function T-i becomes 

(1 + t)c^ 1 log(l - cHx^)dx 



V,{c) = logc / 



7rv4 — 



(1 + ty 1 1 + VI - 4:cH 
= logc -log 

and thus 

The solution to 'H'(c) = with c(0) = 1 is 

c(t) = ^ / . 

2v/2t(l +t) 

From this, using (19), gives (98). 
For (2), we have 



7/(6, c) = logc — — - 



2 4 2 2 y„2 7rV4 



2 



= log(c) log . 

^^^ 2 4 2 2 2 

The critical point (6, c) satisfies the system 



(l-v^b+A/(l-v^'')^-4te2)A/(l"V^^')^-4*c2 

6-V^+ , j} = 0. 

V(l-v^b)2-4te2 



The solution to this system such that c(0) = 1 is given by 



l + 14t + t2-(l+t)Vl-10« + i2 1 _ 6t _ VI - 10« + 42 
= i;^^ = 

Then (19) together with some simplifications give (103). 
For (3) we have 

(1 + t)c^ (1 + t)b'^ bVt 1 log(l - ^/tcx - Vib)dx 



7/(6, c) = logc 



2 4 2 2 7„2 TTv^ 



2 



(l+t)c2 (l+t)62 6v^ 1 1- v^6+ J(l- Vt6)2-4tc2 

= log(c) — ^ ^ log , 

^ 2 4 2 2 2 

The critical point (6, c) satisfies the system 

■ ]^ _ ^2 , 2ic2 _ n 



{l-^/tb+^/(l-^/tb)2-4tc^)^/{l-Vtb)^-4:tc'2 

-(1 +t)b-Vi+ , ^ = 0, 
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with the solution satisfying c(0) = 1, being 



l + 16t + 16t2-(l + 16t)Vl-8t-8i2 i_4t_^i_8t-8t2 

c t = -^^ = and bit) = j= . 

3(l + t)V2t Q{l+t)^/i 

Finally, using (19) one obtains (107). □ 

Next we present algebraic and differential equations satisfied by Qq = J^oit) and the recur- 
sion relation for the Taylor coefficients (/„) of J^o(t) and their exact asymptotic expansions. 

Proposition 8.5. (1) For the potential Vi, Qq satisfies 

~2t + 13t^ + + 2{l-9t + 3f + Abt^ + 32t^) + 64^^(1 + tfg^ = , 
and J-'o(t) satisfies 

(110) t{8 + 7t) - 2 (3 - 4t - 21^2 - Uf) T'^ - 2t (l - 5t - - 7t^) J^o = , 
with J-o(O) = 0, J^o(O) = and (/„) satisfies 

49n'(l + n)fn + 7(1 + nf{32 + 21n)/„+i + {2 + n) (544 + 543n + 139n^) fn+2 

+ {3 + n) (224 + 157n + 33n^) f^+z - 8(2 + n){A + n)(6 + n)/„+4 = 0, 
with /i = 0, /2 = 1/2, = 5/6. For large n, the asymptotics of (/„) is given by 

147 rr 7" / 105 16065 1109115 

(111) tn - ^y^^ V 32^ ^ 2048^2 ~ 65536^3 ^ 

(2) For the potential V2, Go satisfies 

(112) 1 - 13t + 22^2 - 9t^ -t^ + (-2 + 32t - 108t^ + 76t^ + 2t^)go - lOSt^G^ = 
and J-'o(t) satisfies 

3 - 5t + + - 2 (3 - 17t + 5t^ + t^) J'o - 2t (1 - lit + llt^ - t'^) 7^ = 0, 
with Jo(0) = 0, J'^(O) = 1/2 and (/„,) satisfies 

(-2 + n){-l + n)nfn - 2(-l + n){l + n){2 + 5n)/„,+i - (2 + n)(44 + 2bn + bn')fn+2 

+ 4(3 + n) (116 + 91n + 17^^) /„+3 - (4 + r2)(1326 + 701n + mn^)fn+A 

+ 2(5 + n)(6 + n)(85 + 19n)/„+5 - 3(5 + n)(6 + n)(8 + n)/„+6 = 0, 

with fi = 1/2, /2 = 3/4, /s = 8/3, /4 = 12, /s = 312/5. Moreover, for large n, the asymp- 
totics of (fn) is given by 

J_ ,/2 (5 + 2v/6)" / _ 45v^ 8435 _ 238805^6 ^ f J_ 
^ ' ^" 3v/^V3 ^7/2 32n 1024^2 32768^3 ^ 

(3) For the potential V3, Qo satisfies 
(114) 

-2t+llt=^+65t^ + 107tV81t^+27t^ + (2-20t-22t2 + 184t3+560t^ + 700t^+432t^ + 108t^)^o 

+ (108t^ + 540t^ + 1080t^ + 1080t^ + 540t^ + 108t*)^^ = 
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and J-o(t) satisfies 

-8t- 35t^ - Ut^ - 16t^ - 6 (-1 + 15t^ + 34t^ + 28t^ + 8t^) J'q 

-2{-t + 5t^ + 29t^ + 47t^ + 32t^ + 8t^) J'q = 0, 
with J-o(O) = 0, J-'o(O) = and (/„) satisfies 

128^2(2 + n)fn + 32(1 + n) (66 + 97n + 27n^) fn+i + 8(2 + n) (l836 + 1568n + 305n^) fn+2 
+ 4(3 + n) (9972 + 6193n + 929^^) /„+3 + (4 + n) (54276 + 26661n + 3257n2) /„+4 
+ (5 + n) (38220 + 15479n + 1595^2) /„+5 + 3(6 + n) (3972 + 1349n + I2ln'^) /„+6 

+ 5(7 + n) (36 + 5n + n^) /„+7 - 8(6 + n)(8 + n)(10 + n)/„+8 = 

with /i = 0,/2 = 1/2, /s = 3/2, /4 = 47/8, /5 = 138/5, /g = 430/3, A = 11175/14. More- 
over, for large n, the asymptotics of (/„) is given by 



;il5) fn 



(4 + 2v/6)" / 5 (62 - 23^6) 35 (4567- 1858^6) 
' y 8n ^ 64^^2 

35 (2608410 - 1064767^6) 



with 



512/23 

32 



O 




3a/(267+ 109v/6) 

Remark 8.6. Notice here that the exponential rates change depending on the counting prob- 
lem at hand. Excluding just a few types of vertices leads to different exponential behavior. 
It is worth pointing out that in 

J. -n 
fn - C V 

although the exponential growth rate to depends on the details of the model, the exponent 
7 is universal, as was observed in [LGMll] and also in [BBG12, Sec. 2]. 

The recurrence relations for the coefficients (/„) are not in general of the lowest degree. 
However we did not attempt to simplify them even further because they are easily deduced 
from the differential equations satisfied by J-'q. 

Remark 8.7. The linear recursions for (/„) or the linear differential equation for J-q cannot 
compute the Stokes constants, i.e., the leading terms in the asymptotic expansion (111). 

It is the algebricity of /o(^) which uniquely determines the Stokes constants. In the case 
at hand the Stokes constants come from the explicit expressions of J-q. 

Proof. The results from equations (110), (112) and (114) follow straightforwardly from equa- 
tions (98), (103) and (107). 

The idea of proving these asymptotics is based on the analysis of singularities as explained 
in [FS09]. This is particularly tractable as we have explicit expressions for the planar limits. 
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For Vi, from (98) one can see that the singularities of J-q are at to = 1/7 and to = ~1- The 
smallest singularity gives the leading terms. In our case the expansion of J-q near to = 1/''' is 

1, , .^ , 3(l-7t) 137(l-7t)^ 49 [7, „ ,5/2 5569(l-7t)^ 

— -7- 12 bg(2 +8 bg 7 - + ^ \ - 1 - 7tf'^ + ^ - 

16^ &v ;^ B\ jj -r ^ 320 V2^ ^ 12288 



343 



- (1 - 105473 (1 - 7t)' 51793^ | (1 - 7tf' 



+ — .oL^. -7^. + - 7ty). 



1024 131072 98304 

Since the main contribution to the asymptotics of the coefficients is given by the half powers, 
combined with 



(116) (i-x)" = 5^(-i; 



k=0 



applied for a = 5/2, 7/2, 9/2, after a simple asymptotics expansion give the result of (111). 

For V2, the same argument works in this case for the other examples. Using (103), we can 
see that the singularities of J-q are = 5 — 2a/6 and 5 + 2^/6. For the asymptotics of the 
coefficients the leading one is the smallest, namely to = 5 — 2\^. The expansion near to is 
in this case 



12 "'^ 3 K'"^J^^^''"'^H'-w 



' ^ '1-- 9 + 2^6 1-- -— A 1008 + 1270a/- 1-- 



45V3y V to J 27 V J\ to J 63 V v3V to J 



•14 l\f ty y 1089936 + 1337137^|(i-±)'^' /. ^ . 5n 

Considering the coefficients of the half powers, and noting that 1/to = (5 + 2-\/6), one gets 
(113). 

For V3, we proceed similarly. From (107), observe that the singularities of To are to = 
\ (—2 + a/6) and \ (—2 — a/6). The one with smallest absolute value is to = :| (~2 + a/6). 
The expansion near this point is given by 

23 /t: 1. /2/„ ^\\\ 1 /„„ t\ /2539 f2\ ( ^ t^^ 



12 ^+2'°H3('-'^)jj+Is('"-'"'^)i'-W + i-i^-''%3n- to 

/ ,\5/2 / x7/2 

256 1 - 1 . . / . \ 3 1472 1 - ^ 



^ + f 122589 - 50038 Vel ( 1 * * V *o 



45^267 + 109^6 ^4 V ) \ to ) 21\/7929 + 3237\/6 

1 / _x / f\4 16520 (1 -f I // .\5^ 

+ _L ( 5233579 - 2136536^6 ) ( 1 - , ^ ' = + 0(1-- 

216 V /V to) 81V26163 + 10681^/6 VV to, 

from which, noticing that 1/to = 4 + 2a/6 one can deduce (115). □ 
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9. The planar limit for extreme face potentials 
Consider the extreme formal potentials 
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:iir) V7(x) = ^-E 



;ii8) v,(x) = ^-E 



2 ^ 2n 

n.>4 

™2 j.n/2—l^n 



2 ^ n 

n>3 

We consider now the case of extremal potentials and compute the corresponding planar 
limit. 

Remark 9.1. For simplicity, in this section we will drop the subscript / from the writings 
of 7^/, J'o,/, V/. 

Proposition 9.2. (1) For the potential V*^^, the planar limit J^o(t) has the following 10 
terms in the Taylor expansion 

, , t 47^2 49t3 11839t^ 9283^^ 3260543t6 18387797t^ 941448191t^ 
•^o(i) = 77 + -rrr + —r + .on + m + — ^ — + tfz — + 



2 24 4 120 10 336 168 720 

:23647t^ 931715351 
"30 ^ 440" 



490223647^9 93171535189t^° 



and its radius of convergence is to = ^-j^- In addition, if /„ is the coefficient of in the 
Taylor expansion of Fq, then the asymptotic is 
(119) 



„ _1 / 2^-l(i^) A 



243 - 8v^ 91881 - 2640v^- 5696v^ ^/l 



72n 10368n^ V n 



(2) For the potential V, the planar limit J-q has the Taylor expansion 

_ , , 7t 109^2 15631^3 256629t^ 38720767^^ 658811733t^ 

= IT + — + -60- + + ^10— + 112 +^'"- 

The planar limit has radius of convergence given by to, the only positive root of the polyno- 
mial equation 

-11 - 128t + 41088t^ - 20480t^ + 4096t^ = 0. 
The coefficient /„ of from the expansion of has the asymptotic expansion 

(120) /„.e(i/^(l,*^^^0(l 
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with 



/34133 - 914556^0 + 449856^2 - 89344^3 
V 1768687r 

36145645 + 79913928to - 39094848^^ + 7808512^3 
11319552 

7806311269 + 20984001752to - 10129539392^^ + 2006727168t;] 

1026306048 



to = 0.0180827901833 . 
1/to = 55.3012001942... 
c = 0.1786898225... 
di = -3.3197404318 . . . 
d2 = 7.9727292073 .... 



Proof. For part (1), we will find the singularities of 7^(t) = c^{t) and then expand it around 
its singularities nearest to the origin, as explained in [Garll]. 
First, notice that 



and 'H'(c) = implies that the equation satisfied hj TZ = is 



Our condition is that 7^.(0) = 1 and this determines the branch near 0. One can actually 
find the solution explicitly, however that is not very useful. The singularities of TZ are at 
the points where the discriminant of (121) vanishes. That means that the singularities are 
solutions to the polynomial equation 



-16(-5t2 + 96t^ - 96t^ + 32t^) = 0. 



Numerically, 



n{c) = log(c) - - — log 



1 + VI - 4:tC^ 

2 



(121) 



i-n-t^ + mH^ + 47^t^ - wnh"^ + lenh"^ = o. 



The latter are given by 
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The singularity of 7^ is thus to = ^ {2 - ^) ^ 0.0550592. The series expansion of TZ near 
to is 



^^^^ 15 ( 28569 + 23328v^ + 17746v^') ( 1 ^ 



to 



1574640 V V / \ t 



1^0 



_64^ / ^ A _ 1\ 

1574640 V ' \ toj 



2 



1574640 



15 ( 37321489 + 30114648 + 30114648^) f 1 - - 



5/2 



Furthermore, the simplest way to proceed from here is to notice that J-o(t)/t^ differentiated 
three times is TZ'{t)/Tl{t). This in particular means that the singularities of J-q are the same 
as the ones of R and eventually the zeros of TZ. Since the zeros of TZ are only ±1, it follows 
that the singularity of J-q is also to- The expansion of the third derivative of J^o(t)/t'^ near 
to is thus 



1 /I /oo. , „^ 3/t\ , T 424 + 308 + 256 



(^832 + 664 + 1 76 ^) / W 1 



3 V 5 V / V to 135 



/19984+ 16062v^+ 12740v^ f T 12008 + 7336^+ 7472^/ t 

27 V ~ To 3645 V ~t'o 
\/ll3890440 + 92086727v^ + 76015706v^ / t o((l ^ 
1458^2 V to) VV to 

Using (116), we can deduce the behavior of the coefficients of the third derivative of J-o(t)/t^ 
and then a simple exercise leads to the asymptotics of the coefficients of J-q from (119). 
For part (2) we proceed similarly. This time. 



5^ b 1 li-bVt+\/{i-bViy-4tc- 
nib, c) = logic) -c -y-^-^iog 

From the critical system satisfied by b and c eliminate b and then consider 7^ = to arrive 

at the equation satisfied by 71: 

(122) 

1447^4^2 ^7^3t(60- 192t) +7^2 (-2 - 52t + SSt"^) +7^ (l + 15t - IGt"^) - (-1 + 2t - t^) = 0. 
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We are interested here in the branch which at is 1. The singularity points of TZ are at the 
zeros of the discriminant. These are in our case the roots of 

-11 - 128t + 41088t^ - 20480t^ + 4096t^ 

The solution we are interested in is the only solution to in (0,1). Approximately, to ~ 
0.0180827901 .... The value of TZq = TZito) can be found in terms of to as 

_ 1 856to 1280^2 256t|j 

Using Newton's method described in [FS09, VII 7] one can see that the singularity of TZ is 
of square root near to- To find the expansion, write 



To) 



7^(^) = 7^o + X] «fc (i-^l 

fc=i ^ 



where M is the desired level of approximation. Plug this into Theorem 1.4, expand everything 
near to, match the coefficients and then solve the system thus obtained for a^. In our case 
we can take for simphcity M = 3 and solve for oi, 02, 03. The system in this case is of the 
form 



ajvis + V2i = 



2a2f 13 + ajvu + V22 = 

2a3aifi3 + a2fi3 + 3a^a2fi4 + a\vi5 + 021^22 + a^f23 + ^^'31 = 

2a4aifi3 + 20203^^13 + 3ai{al + 0103)^14 + 4af 02^15 + 03^^22 + 2aia2t'23 + alv24 + 01^32 = 
205011^13 + (a| + 2a2a4)t;i3 + {al + 6010203 + 80104)^14 + Qajalvi^ + 4af 03^15 + a4?;22+ 
+0-2^23 + 2aia3i;23 + 3aja2V24 + afv25 + a2V32 + ajuss = 

where the matrix (fij)j=i,3;j=i,5 with coefficients in Q(to) is given in reduced form by 

17+64^ ^2to UAtl " 

-3595+23184*0 -10944tg+2048ig 23-10903to+5440*g-1088fg i r+ i QO+2 qa+ occ+2 
J5g4 33 ~y ~ iOto + oZIq —O^to — ZOOto 

1769-21424to+12352tg-2048fg -2+631to-320fg+64fg 3-64fg , „ , 1/1/1+2 

17424 33 8 "'"'^''O "'-'^'^''O 

There are two different solutions for oi, a positive and a negative one. The appropriate 
one is the negative one in our situation because R has only non-negative coefficients (see 
[BDFG02] for a proof of this). Once ai is solved, the other coefficients are determined 
automatically in a unique way. Also notice here that ai is a square root of a number in 
Q(to), and that all for even k are in Q(to)! while for k odd are in Q(to)/'2i- 

Now given the expansion of TZ near to, the rest follows as in the previous case. Namely, 
we can find the expansion of TZ' (t) /7Z{t) near to and thus the asymptotics of the coefficients 
for TZ'(t)/TZ{t). In turn, since J-o(t)/t^ differentiated three times is exactly TZ'(t)/TZ(t), the 
proof of (120) follows straightforwardly. 

Worth mentioning is the fact that the constant C from (120) is C = —j^^r^^ which explains 
the square root expression of C, while the other constants di and ^2 are in Q(to). □ 
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Remark 9.3. (1) The expansion in (120) can be improved to 

-(l/*o)" 



1=1 

where C is the one from (120) and the constants dn are actually in vi'-oj 
(2) J-'q" is an algebraic function and this determines the Stokes constants in the previous 
result. However the algebraic equation is very long and this is the reason for not 
including it here. In addition, J-'q also is the solution to some algebraic equation, 
though this is very long either. The differential equation satisfied by J-q implies a 
recurrence relation for the coefficients (/„) which is again very long, thus not included. 

10. Other examples of planar limits 

Among other computations we mention here the case of counting planar diagrams with 
vertices of valences 3 or 4. This corresponds to the case of potentials given by 

2 A 

^ ^ 2 3 4 

for the counting of diagrams with a fixed number of edges. The problem of counting planar 
diagrams with a fixed number of faces corresponds to the potential 

2 d 

^ ^ 2 3 4 

The calculations are very similar to the ones for the extreme potentials in Sections 8 and 9. 
The results are as follows. For Vi, the asymptotics of the coefficients of Fq are given by 

where t = to is fhe closest root to 1/5 of the polynomial equation 

= 6912 - 13824t - 146592t2 - 239488t^ - 2602569t^ - 4300752t^ + 79091888^*^ 
+304167552t^ + 410284704t^ - 1349207040t^ - 7615156224t^° 
-4603041792t^^ + 31506516736^^^^ 

and C, di,d2 G Q(to) are given numerically as 

to = 0.2094195368 . . . 
1/to = 4.7751036758... 
C = 1.4826787729... 
di = -7.2166440681 . . . 
^2 = 37.5616277128.... 

Similarly for the potential V2 we have 



j^7/2 I I j^2 
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where to is the closest root to 0.023 of the polynomial equation 

= -43625 - 614400t + 89812992^^ + 895478272t^ - 3041722368t^ - 11466178560t^ 
+32248627200t^ 
In addition, C, and di G Q(to) are numerically approximated as 

to = 0.02305646139... 
1/to = 43.3717899396... 
C = 0.2023938212. . . 
di = -3.2617202693 . . . 



In both cases one can compute the asymptotics of the planar limit in the form 

;iAo)" 



fn = C- 



n 



7/2 



for any M > 1. 



11. Analyticity of the Planar Limit 



In this section we prove Theorem 1.1 and some consequences. Let us introduce some 
notation. For a given sequence a = {«„}.„>! in one of the spaces ^^^(N), define 

il/n 



a[aL) = sup |a„| 

n>l 



Theorem 11.1. (1) For even potentials 



V(-) = y-E 



2n 



n>l 



In 



if a (a) < y/S, then the planar limit J-'Q^(a) is absolutely convergent as a power series 
in infinitely many variables. In particular J-'q is analytic on -B^J^- 
(2) For the potential 

2 n 



n>l 



if Q!(a) < vl2, then J^o{^) (in absolutely convergent series in infinitely many 
variables. In particular J-q is analytic on B^j^. 

Proof. We can write 



•Fo(a) 



n.=l 



5^ c\ax 

\A has only even blocks / 



^^^^ 

n=l \AI-2n 
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where the inner sum is over partitions of size 2n. Note that cx > 0. Now if |a„| < r"/^, then 

cx\ax\ < r" Cx 

X has or 

and 



X\-2n X\-2n 
X has only even blocks A has only even blocks 



cx\ax\ < r" Y 

X\-2n X\-2n 

Hence, in order to compute the radius of convergence we need to compute the radius of 
convergence of the planar hmit J-q for the case of = r"/^. Similarly, for the radius of 
convergence of J^q^ it suffices to look at the case a2n = r"' and a2n+i = 0. 

Now, in these particular cases, according to Proposition 8.3, the radius of convergence of 
-Fq^ is 1/8 while the one of J-q is 1/12. In fact, it is easy to see that the coefficient /„ of t"" in 
J^Q^it) satisfies /„, < 8" and the coefficient /„, of t"' in J-o(t) satisfies /„, < 12". Consequently, 
we have that 

Y < 8" 

X\-2n 

X has only even blocks 



J] CA < 12". 



X\-2n 



Cx < 8" for any partition A of size 2n with only even blocks 



and 

Therefore, 
and in general 

Cx < 12" for any partition A of size 2n. 
Now, a celebrated Hardy and Ramanujan (1918) result shows that the number of partitions 
of size k is asymptotically ^^^e''"\/^'^/^(l + o(l)) ([And98]). From this it follows easily that 
the series J'q^ converges for any r < 1/8 and J-q^ converges for any r < 1/12 which concludes 
the proof. □ 

Given a power series in the form 

2 n 

n>l 



we set 



il/n 



a(V) = sup |a.„| 

n>l 

It is clear that V is analytic near if and only if a(V) < oo. With this notation we have 
the following corollary which confirms t'Hooft's conjecture. For the following statement we 
denote the planar limit J^Q^{t) and J^o{t) to be the planar limits obtained by replacing a„ by 

t I Ciji' 

Corollary 11.1. If V is even then J-'^^it) has radius of convergence at least 



For arbitrary potentials V, J^o(t) has radius of convergence — j= 



12y/^' 

Both of these bounds are sharp. 
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The fact that these bounds are sharp, follow from Proposition 8.2 or Proposition 8.3. As 
made clear from the examples in Proposition 8.4, for the same a(V), the radius of convergence 
can be larger than the one given in this Corollary. 

Remark 11.2. The analyticity of J-q and J-'q^ in infinitely many variables can be deduced 
also from the perturbation result in Theorem 12.1, though without any estimate on the 
radius of convergence. 

Remark 11.3. The reader might wonder what happens with the planar limit if instead 
of considering the potentials V(x) = ^ — J2n>i consider the potentials V(x) = 

^ — X]n>i (^n^"'- In this case the extreme potentials are given by the case where a„ = t"/^ 
and this is V(x) = x^/2 — 1/(1 — x\/t). It turns out after some analysis that the radius of 
convergence of J^o{t) in this case is given by to? which is the only solution in (0, 1) of the 
polynomial equation 

= -226492416 + 962592768t + 34574598144^2 + 334387408896t^ + 7450906184352t^ 
+21095006644064t^ + 130097822364531t'^ + 55792303752096t^ + 67902575063040t^ 
+19100742451200t^ + 6115295232000t^°. 

Numerically this is approximately to = 0.04955391 .... In this case, the planar limit as a 
function of the coefficients a„ is an analytic function on B^. 

Remark 11.4. It would be interesting to know what happens with the case of sequences 
a„ which are not in Apparently the £}. is important for the well definition of convergent 
geometric series in infinitely many variables. 



12. Perturbation Theory 

The main result of this section is a stability result. It says that given a potential whose 
equilibrium measure is one interval, then, under some non- degeneracy assumptions, any small 
perturbation preserves the one interval support of the equilibrium measure and in addition, 
the planar limit depends nicely on the perturbation. 

Before we state the result, we want to define a class of perturbations of a given potential 
V. This definition is long and depends on many parameters, however the idea is quite 
simple. We want to take perturbations of V so that the maximizer of the function F can 
be parametrized in a nice way. The reasonable way of doing this is to have perturbations 
close to V on some open interval containing the support of fiy and large outside this open 
interval. 

For this purpose, assume that X is a Banach space over the reals which will be the ambient 
space of the parametrization. Now, given an open subset D of X such that e D, and /, J 
open sets of M, an integer k > 1 and R,6>0, we define U{k,V,D, I, J, R,S) the class of 
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functions V : D x M — )■ M in two variables, with the properties, 

1) V{0,x) = V{x) 

2) for each t G D, x — V(t, x) satisfies (5) 

3) (t, x) V(t, x) is C^'3(D X M) 



(123) 



4) sup \V{t,x) -V{x)\ < 6, k sup ||Hess^(V(t, ■) - < ^, 

t6D,a;e/ t6D,2;eJ 

5) inf (V(t,x) -21og|x|) > i? 

t6D,a;^/ 



where, C^'^ stands for the set of jointly differentiable functions in (t,x) with k continuous 
(Frechet) differentials in t and three continuous derivatives in x. Also, HesSa; stands for the 
Hessian with respect to the variable x and || ■ \\hs is the Hilbert-Schmidt norm. 

In words, 1), 2) and 3) of (123) define the perturbation which is assumed of class C*''^, 
while 4) means that the perturbation is uniformly close to ^ on D x / while the Hessians are 
uniformly closed on D x J and 5) encodes the fact that outside the interval /, the perturbation 
(minus the logarithmic term) is larger than a constant R uniformly in the parameter t G D. 
We introduce here the interval J because as we will see below in the proof of Theorem 12.1, 
we only need the Hessians close for x on a neighborhood of the support of 

The reason of introducing condition 5) in (123) instead of condition 4) with / = M is 
because for large values of x, we do not need the perturbation to be close to V . We only 
need the perturbation to be large for large x. Actually, 4) and 5) constitute a weakening of 
the condition that the perturbation stays close to V uniformly on the whole M. 

Recall that we set 

fe>W^r (^'(" + '')-'^;(f^''))''^ v.. [-2,2], 

J-2 {x - y)'K A - y'^ 

Theorem 12.1. Assume that V is aC^ potential satisfying (5) with H{c, h) and ipc,b defined 
by (57) and (56) respectively. Suppose that the following conditions hold true: 

1. {c,b) is the unique maximizer of H; 
"^"^^^^ 2. ipc,b{x) > 0, for all x E [-2,2]. 



Under these assumptions, there exist 

• an interval / C M, 

• positive numbers Rq and Sq 

with the property that for any choice of 

• R > Ro, < 6 < 6o 

• an open neighborhood J of [—2c + b,2c + b], 

• a Banach space X 

• andV eU{k,V,B,I,J,R,S), 
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the following hold 

1) there exists an open Dq C D with G Dq and 

2) (c, 6) : Do (0, oo) x M which is such that c(0) = c, b{0) = b 

3) (c(t),6(t)) is the unique maximizer of H{t, ■) ( defined by (57) for V(t, ■)) 

4) Do X [-2,2] 3 (t,x) '?/'c(t),fe{t)(a;) G M zs positive. 

Furthermore, the equilibrium measure for V{t, ■) has a single interval support for t G Dq 
and the planar limit F^^t = -fo,v(t) is a function on Dq. 

In addition, if X is either a finite dimensional space or of the form X = {(an)n>i C M : 
X]n>i knk" < oo} for some r > and V is real analytic on a neighborhood of the support 
of fly such that (t,x) — V(t,x) is real analytic on a neighborhood of x [— 2c + b,2c + b], 
then, we can take Dq so that t — c(t) , t — )■ b{t) and t — )■ F^^t are real analytic functions. 

Proof. The key point of the proof is the fact that the maximizer (c, b) of H is unique and 
isolated and then by perturbing a httle bit the potential V, the maximizer of H{t,-) is 
to be found near (c, 6). Finding the maximizer (c(t),6(t)) of H{t,-) boils down to finding 
the critical point of this function near (c, 6). This can be achieved by the implicit function 
theorem and the fact that the Hessian of H is non-degenerate near (c, b). 

Now technicalities. The first thing we want to do is to prove that for the unperturbed 
function H, (c, b) is a non-degenerate critical point. To do this we want to check that the 
Hessian of H at (c, b) is positive definite. For simplicity of the discussion, we will assume 
without any loss of generality that c = 1 and b = 0. Now the non-degeneracy is equivalent 
to the fact that 



(126) 



r2 x^V"{x)dx r2 xV"{x)dx 

J-2 ttV^^ J-2 -Ks/t^ 
1-2 xV"{x)dx 1-2 V"{x)dx 

J— 2 nV^—x'^ J —2 7rV4— a;2 



is positive definite. 

Recall that the critical point equations give 



2 xV'{x)dx ^ , V'{x)dx 
— 2 and ' 



_2 7rV4 — x^ 7-2 nyj A — x"^ 



Integrating by parts the first of these one deduces that 

^ x'^V"{x)dx ^ V"{x)dx 



2 n^/A — J^2 TT\^A — x"^ 



Armed with this, the non-degeneracy of the Hessian (126) follows once we prove the following 



(127) / ^ — >0. 



2 



[2 ± x)V"{x)dx 



T\\/ A — x^ 



This follows from 

[2±x)V"{x)dx C <i nn, X ^ 



/ ^{v'{x)-v\±2)y 

J-2 dx 

2 {V'{±2) - V{x))dx 
-2 (±2 - x)7rV4 - x2 



il^{±2) > 0. 
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Let M = H{c,b) be the maximum of H. For any choice of e,r > with r > e > 0, 
obviously one has 

(128) sup{H{u, v) : u > 0,v e R,r^ > [u - cf + {v - hf > e^} < M. 

Indeed if this is not the case, then there is a sequence (c„, 6„) such that > (n„, — c)^ + 
{vn — by > so that limn^oo H{un,Vn) = M. Passing eventually on subsequences, we may 
assume that n„ and Vn converge to u and v. Clearly u ^ 0, otherwise H{u, v) = —oo. This 
implies that H{u, v) = M and at the same time, (m, v) is within positive distance from (c, 6), 
hence contradicting the uniqueness of the maximizer. 
Next, consider 

U{x) := V{x) - 21og|x|, 

and notice that from (5) we have U{x) > C > — oo and lim|a;|^oo U{x) = +oo. Now we use 
(37) to justify that 

/"^ U(ux + v)dx \og\x + V /u\ Uiux + v)dx 

H{u,v) = - / , / , < - 



27rv/4 - x^ J-2 27rV4 - x^ ~ J-2 - x^ 

Assuming that (u — c)^ + (t' — 6)^ > r^, it is easy to deduce that, > ^Jr"^ jl — — 6^ 

and thus, 



(129) H[u,v) < 



^ U{ux^v)dx f'^U{ux + v)dx f ^U{ux + v)dx 



27rV4^x2 Ji 27rV4 - x"^ J-2 2ti^/A - x'^ 



< -C/3-h{^/ry2-c^ -b^) 

where h{x) = in{\y\>x U{y)/6. In particular, for large r we learn that H{u,v) < M. 

Equations (128) and (129) guarantee that for any e > 0, there exists 6o G (0, 1) such that 

(130) H{u, v) <M- 35o 

for all {u,v) outside a ball of radius e around (c, 6). We take Rq > such that |x| > Rq 
implies h{x) > -C/3 - M + 3 and define / = [-Rq, Rq] U [-2c - 1 + 6, 2c + 1 + 6]. The 
purpose of this choice of / is to make it a neighborhood of the support of fiy 

With these choices, for any R > Rq, < 6 < 6q and V G U{k,V,'D, I, J, R,6), from 
the conditions 4) and 5) of (123), and the reasoning which led to (129), one gets for r = 
^J2{RI + 62 + c2) that 

1. \H{t,u,v)- H{u,v)\<5, for r"^ > {u - cf + {v - bf > 

2. H(t, u,v) < M — 3 for (n — c)^ + {v — by > r^. 

We are led to the conclusion that for all t G D, max„>o,t,ei8{-f:r(t, u, v)} is attained for (m, v) 
in the ball of radius e around (c, 6). Indeed, otherwise, assume that there is a maximizer 
(m, v) outside the ball B^{c, b). Since, \H{t, u, v) —H{u, v)\ < 6, combined with (130), implies 
that H(t,u,v) < M — 25. This contradicts 4) of equation (123) from which we gather that 
if(t, c, b) — H{c, b) > —6, or H(t, c,b) > M — 6 > H{t, u, v) + 6, thus (m, v) can not be a 
maximizer of H{t,-). 

The maximizer is a critical point of H{t,u,v), therefore Vu.vH{t, ■) = 0. To solve for 
{u,v), we interpret it as the definition of an implicit function t — t- (c(t),fe(t)). This can be 
done thanks to the combination of the last part of 4) of (123), 1) of (124) and the implicit 
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function theorem. These yield for a set Dq C D, which contains that there exists a C'^ 
function t — > (c(t),fe(t)) which is the maximizer of H{t, ■). Taking a smaller subset of Dq, 
it is easy to show that 'ipc{t),b[t) > on [—2, 2] and the dependence of Fo^t on t is a simple 
consequence of (51). 

In the case of analytic perturbations with X a finite dimensional space, the only thing 
we need to point out is that (cf. [KP02]) the implicit function theorem produces analytic 
versions c(t) and b{t) for t in an eventually smaller Dq. The analyticity of F^t follows from 
(51). 

On the other hand in the case X = {(a„)„>i C M : J2n>i I'^nk" < oo}, one needs a bit 
more work. The analyticity of functions in infinitely many variables is trickier than the case 
of analytic functions in finitely many variables. However, our space here is essentially 
over the real numbers and for this case many things are like in the finite dimensional cases. 

What we mean here is that for the case of i^^N) over the complex numbers, the theory 
of analytic functions is treated in [Lem99] and [Rya87]. The main results are that every 
holomorphic function on £^(N) has a power series expansion and every absolutely convergent 
power series expansion defines a holomorphic function. 

In our situation, the functions are real analytic (meaning they have a power series ex- 
pansion), thus by complexification they become complex analytic and therefore they are 
holomorphic functions. Then, for the complexification, we know that the implicit function 
theorem yields that the resulting functions c(t), b{t) and Fq i all are smooth functions of t 
on a small neighborhood of X. Furthermore, since F is actually a holomorphic function it 
is not hard to prove that the choices of c(t), b{t) and Fo^t can be made holomorphic. Using 
this we can conclude that the real parts of c(t), 6(t) and Fo,t are real analytic. □ 

Appendix A. The first few terms of TZ, S and J'q 

In this appendix we give the first few terms of the unique solution (i?, S*) G ^ of Equations 
(13) and also of the formal planar limit J-q. 

Q O O O 

S = ai + aia2 + 203 + 0102 + a^^oa + 40203 + 601O4 + 605 + 01O2 + 30^^0203 + 602O3 + 801O3 

o 4. 

+a]^04 + 18010204 + I80304 + 120^05 + I802O5 + SOoiOe + 2O07 + 01O2 + 60^0203 + 802O3 

"-iO 0'?'? 00 Ozl 

+20;^03 + 32O1O2O3 + I2O3 + 40]^0204 + 3601O2O4 + 420^0304 + 72O2O3O4 + 54O1O4 + 0^05 

+4801O2O5 + 3602O5 + 108010305 + 7204O5 + 20of 06 + I2O01O2O6 + 8O03O6 + 900^07 

+800207 + 140O1O8 + 7009 + 0(0^^) 

TZ = 1 + 02 + 02 + 201O3 + 3o4 + 02 + 601O2O3 + 4o| + 3o^04 + 902O4 + I201O5 + lOog + 02 

22222 2 23 

+12010203 + 601O3 + 160203 + 12010204 + 180204 + 42010304 + 1804 + 40105 + 48010205 

O ^QOOOOQ 

+360305 + 30O]^O6 + 4O02O6 + 6O01O7 + 35o8 + 02 + 2O01O2O3 + 30o;^O2O3 + 4O02O3 + 3201O3 

fjfjQo 0000*? 

+30o;^O2O4 + 3O02O4 + 200^0304 + 2IO01O2O3O4 + 8403O4 + 63o;^04 + 9O02O4 + 2O0J02O5 

A 

+I2O01O2O5 + I32O1O3O5 + I8O02O3O5 + 252O1O4O5 + 72O5 + 50i06 + I5OO1O2O6 + IOOO2O6 

+26O01O3O6 + I5O04O6 + 60ofo7 + 3OO01O2O7 + I6O03O7 + 210o^O8 + I7502O8 + 28O01O9 

+126010 + 0(0") 
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Ci 02 1 2 

y + y + 2«i«2 



^2,l22|^2|l23i^2,l24|^2| i^3 

— Oi an H — Oi H — ai a, H h- ai H — fli 



1 



+2010203 + 0^0203 + 301O2O3 + 20^0203 + 401O2O3 H — + 201O3 + -0^03 + 20203 



+801O2O3 + 4o2o| + 401O3 + y + + + 0204 + ^010204 + O1O2O4 + ^0304 



,3 , "4 



1 



9 



2 2 3 f-T 3 ^ 2 iyti^ 2 2 2 

+901O2O4 + 202O4 + 6010304 + 70^0304 + 2401020304 + 603O4 H — - — h 90^04 + —0204 

2 

1 



-2oi05 + 2afo5 H — a?a5 + 601O2O5 + 8o?0205 + 12010205 + 80305 + I801O3O5 + I202O3O5 
5 



-I801O4O5 + 



I8O5 5o6 



502O6 



^ + — +50i06 + -O1O6 + „ 

2 2 



+ 2O01O2O6 + 502O6 + 2O01O3O6 



T Qj 35 

-604O6 + 5O1O7 + 10ofo7 + 20aiO2O7 + 803O7 H ^ h y OiOg + 70208 + 14oi09 



21010 



+ 0(0 



11^ 



where each is given the degree k. For example the monomial 010304 has degree 2x1 + 
2 X 2 + 4 = 10. The meaning of 0{a^^) is that the degree of the remaining terms is at least 
11. 
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